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ABSTRACT1’ 

Cratering  programs  and  data  resulting  from  numerous  single-charge 
explosion  tests  are  summarised  and  compiled  in  tabular  form.  Analyses 
are  performed  on  these  data  to  provide  means  of  predicting  basic 
cratering  parameters.  Prediction  equations  were  developed  by  use  of 
the  method  of  least  squares.  Means  of  updating  these  tabulations  and 
analyses  on  a regular  basis  by  automatic  data  processing  are  discussed. 

Data  are  grouped  so  as  to  account  for  the  factors  which  primarily 
affect  crater  size  and  shape:  y*eld,  burst  geometry,  and  cratered 

medium.  The  influence  of  other  conditions,  such  as  soil  moisture, 
layered  media,  etc.,  is  also  considered.  Qnphasis  is  on  single-charge, 
dry-land  experiments,  vhich  best  permit  isolation  of  the  factors  con- 
tributing to  the  basic  parameters.  Howeier,  effects  of  environmental 
influencer.,  unusual  charge  geometries,  and  other  factors  significantly 
affecting  craters  are  also  briefly  considered.  Similarly,  ba3lc  ejecta 
information  is  included. 

Trends  in  crater  dimensions  are  shown  by  means  of  graphs  normalized 
to  charge  sizes  commensurable  to  large  chemical  and  small  nuclear  yields. 
Seeding  as  a prediction  tool  is  discussed. 
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PREFACE 


This  report  deal3  with  cratering  from  explosive  charges.  It  is  in- 
tended to  compile  in  tabular  form  all  single-charge  cratering  data  suit- 
able for  analysis,  taken  from  numerous  test  programs  conducted  in  a vari- 
ety of  media,  and  it  includes  a somewhat  abbreviated  empirical  analysis 
of  these  data.  The  study  was  conducted  for  the  Office,  Chief  of  Engi- 
neers, Department  of  the  Arny,  under  Task  04,  Project  No.  1»A062118a880, 
by  the  U.  S.  Arny  Engineer  Waterways  Experiment  Station  (WES)  during  the 
period  January  1970  through  April  1971.  The  research  was  under  the  gen- 
eral supervision  of  Mr.  G.  L.  Arbuthnot,  Jr.,  Chief,  Weapons  Effects  Lab- 
oratory, and  Mr.  J.  N.  Strange,  Chief,  Engineering  Research  Branch,  and 
under  the  direct  supervision  of  Mr.  A.  D.  Rooke,  Jr.,  Chief,  Earth  Kinet- 
ics Section.  This  report  was  prepared  by  Messrs.  Rooke,  B.  L.  Carnes, 
and  L.  K.  Davis.  Assistance  in  the  search  for  data  and  in  the  biblio- 
graphical compilation  was  given  by  Mr.  J.  A.  Conway;  Mr.  S.  B.  Price  as- 
sisted in  preparation  of  the  graphs.  Assemblage,  typing,  and  proofing 
of  the  report  draft  were  by  Miss  Virginia  Mason  and  Mrs.  Dean  McAlpin. 
Additionally,  special  assistance  in  certain  subject  areas  va3  provided 
as  follows:  crater  ejecta  - Mr.  J.  W.  Meyer;  charge  stemming  - 

SP5  K.  L.  Knudson;  and  crater-cavity  formation  - 1LT  H.  D.  Hardcastle. 

A previous  WEE  compendium  on  cratering  was  used  as  the  primary 
source  of  cratering  data  prior  ,o  i960,  and  a previous  analysis  of 
crater  data  served  as  a guidelin#  for  this  analysis.  The  usefulness  of 
these  research  efforts  in  preparing  this  report  is  gratefully 
acknowledged. 

COL  Levi  A.  Brown,  CE,  BG  Ernest  D.  Peixotto,  CE,  and  COL  G.  H. 

Hilt,  CE,  were  Directors  of  the  WES  during  the  preparation  of  this  re- 
port, and  Mr.  F.  R.  Brown  was  Technical  Director. 
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NOTATIONS,  ABBREVIATIONS,  AND  DEFINITIONS 


Notations 

C An  observed  cavity  proportionality  constant 

d Crater  depth  (general) 

dft  Depth  of  apparent  crater  at  ground  zero 

Depth  to  limit  of  plastic  deformation  at  ground  zero 

dp  Depth  to  limit  of  rupture  zone  at  ground  zero 

d^  Depth  of  true  crater  at  ground  zero 

c A factor  which  represents  the  enhancement  of  crater  dimensions 
in  multiple  explosions 

h Height  of  apparent  crate.'  lip 

k A constant  used  in  scaling  relations 

K A constant  representing  the  fraction  of  ejecta  volume  contribut- 
ing to  the  formation  of  the  crater  lip 

L Hie  length  of  a row  crater 

n Crater  scaling  exponent 

N The  number  of  charges  in  a row  shot.  Shot  number  in  successive 
charges  fired  on  a common  vertical  axis 

r Crater  radius  (general) 

r^  Radius  of  apparent  crater 

rc  Radius  of  true  crater  cavity 

r ^ Radius  of  charge 

r^  Radius  to  outer  limit  of  plastic  deformation 

rft  Radius  to  outer  limit  of  ejecta 

r^  Radius  to  maximum  lip  height 
r^  Radius  to  outer  limit  of  apparent  lip 

rp  Radius  to  outer  limit  of  apparent  rupture 

r^  Radius  of  true  crater 

Radius  to  point  of  maximum  upthrust 
Range  (distance)  from  ground  zero 
s Charge  spacing  in  a row  shot 

u Height  of  upthrust 

v Crater  volume  (general) 

vft  Volume  of  apparent  crater 
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vc  Volume  of  crater  due  to  compression  (compaction) 

Preshot  volume  of  material  dissociated  by  the  explosion 
ve  Volume  of  ejecta 

vcxp  Volume  of  joint  expansion 

Vj,  Volume  of  crater  due  to  plastic  flowagc  of  the  medium 
v_  Volume  of  fallback 

ID 

vt  Volume  of  crater  lip 

v Volume  of  true  crater 

vu  Volume  of  upthrust 

Vf  Void  ratio 

W Charge  weight 

Z Scaled  depth  o:-  height  of  burst  (negative  if  DOB) 

or  An  exponent  based  on  '.he  adiabati.:  (^ 3)  expansion  coefficient 

7 Unit  weight  of  the  cratered  medium 

4 Areal  density  of  ejecta 

v Dyramic  viscosity  of  the  cratered  mediuM 

p Density  of  the  cratered  u.  ••Hum 

o Compres , shear,  and  r>n-.\e  strengths  or  elastic  properties 
of  the  era''-'*'.!  .aedium 

^bb-ovjatlons 

AN  Amaomur.  nit  ; vj 

ANPO  Amnsorii.  « oit.*ofe/fuH  oi) 

DOB  Depth  of  burst  Ho  center  of  gravity  of  charge)  below  original 
ground 

GZ  Ground  zero,  the  hypocenter  or  epicenter  of  the  burst 
HE  High  explosive 

HOB  Height  of  burst  (to  center  of  gravity  of  charge)  above  original 
ground 

kt  Kilotcn 

Mt  Megaton 

NE  Nuclear  explosive 

NM  Nitrone theme 

TNT  Trinitrotoluene 
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'Common  dimensional  abbreviations  are  used  in  accordance  with  "Wea- 
pon Test  Reports  Preparation  Manual";  DASA-26,  September  1966;  Defense 
Atomic  Support  Agency,  Washington,  D.  C. 

Additional  abbreviations  are  used  in  Chapters  1 and  2 to  designate 
agencies,  and  in  Append!.:  A to  identify  explosion-effects  tests  and 
programs.  All  such  abbreviations  are  identified  where  used. 


Definitions 


Apparent  crater 
Crater  lip 
Ejecta 
Fallback 

Multiple  explosion 
Nail  driving 

Optimum  DOB 
Row  crater 
Row  shot 
Rupture  zone 

Stemming 

Surface  tangent 
(above  and  below) 


The  visible  crater,  bounded  at  the  top  by  the 
original  ground  surface  elevation 

The  region  of  continuous  ejecta  surrounding  a 
crater 

Earth  material  permanently  ejected  from  the  crater 
void  by  the  explosion 

Material,  dissociated  by  the  explosion,  which  has 
fallen  back  within  the  true  crater  void 

The  detonation  of  two  or  more  charges  with  suf- 
ficient simultaneity  and  proximity  to  cause  inter- 
action in  crater  formation 

A blasting  technique  using  successive  explosions 
on  a vertical  axis,  with  each  charge  being  em- 
placed in  the  center  of  the  crater  of  the  pre- 
ceding shot 

The  depth  of  burst  at  which  the  largest  desired 
crater  dimension  occurs 

A crater  or  channel  formed  by  the  detonation  of 
charges  emplaced  in  a row -shot  geometry 

A multiple  explosion  with  the  charges  emplaced 
in  a linear  array 

Material  below  and  beyond  the  true  crater  which 
has  sustained  significant  physical  damage  (frac- 
turing, crushing,  shearing,  etc.)  as  a result  of 
the  explosion 

(verb)  The  backfilling  of  the  charge  emplace- 
ment hole  of  an  underground  charge;  (noun)  back- 
fill material 

A charge  geometry  with  the  surface  of  the  spheri- 
cal charge  tangent  to  the  ground  surface.  (Above 
indicates  the  charge  is  resting  on  the  ground  and 
below  indicates  the  charge  is  buried  one  charge 
radius) 
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True  crater 


The  boundary  of  the  crater  representing  the  limit 
of  dissociation  of  the  medium  by  the  explosion 
(the  crater  prior  to  fallback) 

True  surface  burst  A charge  geometry  with  the  center  of  gravity  of 

the  charge  at  the  ground  surface 

Upthrust  Material  that  has  been  permanently  displaced  above 

the  original  ground  surface,  but  not  dissociated 


CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 


British  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  units  as  follows: 


Multiply 

By 

To  Obtain 

feet 

0.30U8 

meters 

square  feet 

0.09290 

square  meters 

cubic  feet 

0.02882 

cubic  meters 

pounds 

0.1*536 

kilograms 

short  tons 

0.9072 

metric  tons 

pounds  per  cubic  foot 

16.02 

kilograms  per  cubic  meter 
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CHAPTER  1 


INTRODUCTION 


1.1  BACKGROUND 

The  importance  of  explosive  cratering  was  greatly  increased  when, 
through  the  advent  of  nuclear  weapons,  it  became  possible  to  release  in- 
stantaneously enormous  amounts  of  energy  from  essentially  a point  source. 
Certainly  the  best  method  to  establish  the  effects  of  nuclear  explosives 
(NE)  is  through  full-scale  testing.  However,  because  of  the  hazards  in- 
volved and  the  restrictions  imposed  by  international  treaties,  not  to 
mention  costs,  full-scale  testing  is  presently  impractical.  Often,  the 
best  method  to  approximate  nuclear  explosives  or  weapons  effects  is 
through  chemical  or  high-'  ^ *7”)  testing,  despite  difficulties  in 
scaling  HE  results  to  NE  charge  ylcl'- . 

Scaling  difficult  ■ , which  have  absorbed  much  of  the  research  ef- 
fort in  this  field,  staa  j,  .‘imarily  from  a lack  of  definitive  information 
regarding  the  relative  importance  of  various  crater-forming  mechanisms 
and  the  physical  factors  bearing  on  these  mechanisms,  as  well  as  from 
similarity  violations  which  inevitably  occur  in  this  type  of  experimenta- 
tion. They  are  compounded  by  basic  differences  in  energy  release  and 
partitioning  between  HE  and  NE  charges.  However,  by  considering  these 
facts  and  weighing  the  scaled  HE  data  against  the  actual  NE  data  that  are 
available,  the  former  can  be  useful  in  predicting  the  effects  of  nuclear 
devices  or  weapons,  particularly  when  the  nuclear  yields  in  question  are 
below  the  megaton  level. 

Future  weapons-effects  research  will  undoubtedly  involve  the  appli- 
cation of  near-surface  or  belov-surface  HE  and  NE  detonations.  A knowl- 
edge of  the  effects  of  cratering  and  its  associated  debris-ejection  mech- 
anism will  be  a prerequisite  to  any  such  application.  The  crater,  its 
surrounding  zones  of  subsurface  deformation,  and  the  ejecta  field  repre- 
sent varying  degrees  of  damage  or  possible  damage  to  structures  located 
therein  (Figures  1.1  and  1.2).  Moreover,  the  size  of  the  crater  produced 
in  such  an  explosion  is  an  indicator  of  the  total  energy  coupled  into 
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the  medium;  therefore,  the  prediction  of  crater  parameters  is  of  major 
concern  in  weapons  employment.  Crater  size  has  also  been  shown  to  be 
useful  in  normalizing  other  explosion  effects  phenomena,  e.g.,  ground 
motion. 

The  military  applications  of  cratering  research  include  the  capa- 
bility of  weapons  to  damage  or  destroy  hardened  defense  installations, 
to  create  obstacles  or  barriers  in  various  situations,  or  to  provide 
expedient  means  of  excavation.  The  civil  applications  are  mostly  in- 
volved with  excavations  of  canals,  harbors,  etc.,  river  diversion  or 
damming,  underground  stimulation  of  mineral  production,  underground 
storage,  etc. 

1.2  PURPOSE 

In  the  past,  efforts  to  analyze  and  correlate  cratering  data  have 
met  with  considerable  difficulty  because  of  the  large  number  of  reports 
in  which  the  data  are  presented,  and  because  of  the  fact  that  cratering 
data  are  often  secondary  to  the  main  purpose  of  the  research.  The  com- 
pendium of  i960  (Reference  1)  and  enalysis  of  1961  (Reference  2)  allevi- 
ated this  situation;  however,  much  cratering  research  has  been  done 
since  that  time,  and  a fresh  look  at  the  problem  of  tabulation,  correla- 
tion, and  analysis  is  in  order. 

The  purpose  of  this  report  is  to  compile  and  analyze  all  available, 
useful  cratering  data  (both  HE  and  NE)  in  one  report,  and  to  present  it 
in  such  a manner  that  it  will  serve  as  a guide  both  for  cratering  appli- 
cations and  for  planning  future  cratering  research.  It  is  intended  that 
the  compilation  be  in  such  form  as  to  permit  continuous  updating  as  ad- 
ditional cratering  experiments  are  performed.  The  data  are  grouped  to 
permit  ijelation  and  quantification  of  factors  which  significantly  in- 
fluence crater  size  and  shape.  These  include  charge  yield,  shot  geom- 
etry, and  properties  of  the  cratered  medium,  as  well  as  other  factors 
that  are  less  dominant  in  their  influence  on.  the  cratering  process.  The 
purpoce  of  the  analysis  is  to  provide  prediction  techniques  so  that  all 
important  cratering  phenomena. can  be  predicted  with  a reasonable  degree 
of  accuracy. 
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It  is  also  the  purpose  of  this  compendium  to  present  all  important 
aspects  of  cratering  phenomenology.  True  crater  dimensions  (see  Fig- 
ure 1.2)  are  important  as  the  limit  of  dissociation  of  the  medium;  it 
is  unlikely  that  any  structure  can  survive  within  this  boundary.  Both 
true  and  apparent  crater  volumes  are  important  from  the  standpoint  of 
excavation,  and  apparent  crater  lip  height  and  the  interior  angle  of 
slope  are  considerations  in  the  creation  of  crater  barriers  and  the  slope 
stability  of  the  crater  walls.  Also  important  in  the  design  of  hardened 
structures  are  the  limits  of  subsurface  deformation  surrounding  a crater. 
Few  observations  are  available  on  these  limits,  but  they  are  summarized 
in  Chapter  3* 

1.3  SCOPE 

Although  this  report  is  intended  to  describe  cratering  phenom- 
ena in  general,  it  has  been  necessary  to  limit  the  contents  in  some 
cases.  Thus,  single-charge  craterB  produced  in  conditions  of  uncompli- 
cated media  and  on  more  or  less  level  topography  are  considered  at 
length,  while  multiple-charge  arrays  and  craters  occurring  under  less 
usual  geometric  or  environmental  conditions  are  examined  only  briefly. 

It  has  also  been  necessary  to  restrict  the  lower  limit  of  charge  size 
to  1 pound. * Successful  experiments  with  smaller  charges  have  been  con- 
ducted, but  too  often  cratering  results  are  open  to  question,  due  pri- 
marily to  difficulties  in  scaling  medium  properties  so  as  to  be  suitable 
to  such  small  yields.  Data  from  shaped  charges  have  also  been  omitted, 
since  their  use  is  highly  specialized.  Cratering  data  in  certain  un- 
usual media  or  conditions  (e.g.  snow  and  ice,  earth-water  interfaces) 
are  considered  only  briefly  in  this  report,  but  bibliographical  refer- 
ences for  these  experiments  are  included.  The  same  is  true  for  craters 
formed  by  unstemmed  or  partially  stemmed  charges.  In  such  cases,  pre- 
vious U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  reports  are 


* A table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  is  presented  on  page  12. 
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cited  which  provide  compendia  of  these  data. 

A further  limitation  in  content  results  from  a self-imposed  require- 
ment that  this  publication  be  unclassified,  making  it  available  to  the 
maximum  number  of  users.  However,  while  classified  data  are  excluded 
per  sc.  they  have  been  included  in  the  construction  of  the  cratering 
curves  and  the  empirical  expressions  derived  therefrom. 

In  order  to  bring  the  compendium  to  a conclusion,  it  was  necessary 
to  establish  a cutoff  date  for  consideration  of  new  test  results.  This 
was  set  as  the  end  of  calendar  year  1970.  Important  tests  have  taken 
place  since  that  time,  and  these  vill  hopefully  t«  Included  in  future 
updating  of  this  report. 


I 
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(Courtesy  of  Lawrence  Radial  .on  Laboratory 
and  U.  S.  Atomic  Energy  Commission ) 

Figure  1.1  Aerial  view  of  a typical  crater  formed  by  a low-yield 
nuclear  device  at  near-optimum  depth  of  burial  in  basalt  (Event 
Danny  Boy).  Note  the  size  of  boulders  in  the  crater  and  crater 
lip  as  compared  to  the  vehicle  (arrow). 


18 


Figure  1.2  Typical  half-crater  profile  and  nomenclature  for  a buried  charge.  Profiles  and  dimen- 
sions are  symmetrical  about  the  centerline*  Various  radial  and  depth  dimensions  are  indicated  by 
r and  d , respectively.  For  explanation  of  terms , see  pages  9 through  1’ . A crater  formed  by  a 
near-surface  explosion  is  similar  in  profile,  but  usually  without  a distinct  cavity  radius  rfi  . 
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CHAPTER  2 

SYNOPSIS  OF  CRATERING  RESEARCH  EFFORTS 

2.1  HE  TEST  PROGRAMS:  PRIOR  TO  1962 

The  first  significant  efforts  to  investigate  the  phenomena  of  cra- 
tering  from  explosives  appear  to  have  originated  during  World  War  II. 
Both  the  U.  S.  War  Department  and  the  Ministry  of  Home  Security  in 
England  conducted  investigations  in  which  various  types  of  American, 
British,  and  German  bombs  were  statically  detonated  in  different  types 
of  soil.  An  early  study  involving  the  modeling  of  large-scale  cratering 
effects  by  using  small-scale  charges  was  conducted  by  the  U.  S.  Army 
Ballistic  Research  Laboratories  to  estimate  the  amount  of  material  re- 
quired to  fill  bomb  craters.  References  for  these  and  other  tests  sum- 
marized in  this  chapter  are  contained  in  the  Bibliography  (Appendix  C). 

Near  the  end  of  World  War  II,  the  Underground  Explosion  Test  (UET) 
program  was  conceived  to  study  design  criteria  for  underground  struc- 
tures to  resist  the  effects  of  underground  explosions.  The  size  of  the 

largest  test  charge  detonated  in  the  program,  320,000  pounds,*  indicates 

• 

that  the  possibility  of  attack  by  a nuclear  weapon  was  now  also  con- 
sidered as  a threat.  The  UET  program  was  conducted  in  two  major  phases: 
a preliminary  program  of  detonations  in  rock  by  the  Colorado  School  of 
Mines  in  19W5-49,  and  a more  extensive  program  by  Engineering  Research 
Associates,  Inc.,  from  19^9  through  1951.  Although  several  series  of 
tests  were  fired  for  crater  investigations  alone,  the  majority  of  the 
tests  involved  target  structures  in  earth  and  tunnels  in  rock. 

A second  major  test  program  initiated  by  the  postwar  interest  in 
protection  against  bombing  attacks  was  the  Isthmian  Canal  Studies  (ICS) 
conducted  in  19^7 • This  program  vas  conducted  by  the  Canal  Zone  Divi- 
sion of  the  Army  Corps  of  Engineers,  and  involved  the  study  of  cratering 


* Unless  otherwise  specified,  charge  weights  are  given  in  terms  of  tri- 
nitrotoluene (TNT)  or  its  energy  equivalent. 
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damage  and  the  initiation  of  slope  failures  from  bomb  explosions  in  the 
embankments  adjacent  to  the  Panama  Canal.  TOT  charges  ranging  from  6 to 
200  pounds  were  fired  in  five  different  media  representative  of  the  soil 
and  rock  types  along  the  canal  route. 

The  first  nuclear  test  series  at  the  Nevada  Test  Site  (NTS)  in  which 
cratering  was  studied  per  se  was  Operation  Jangle.  The  operation  re- 
quired a limited  HE  lest  series  to  permit  some  reasonable  basis  for  pre- 
dicting the  nuclear  test  phenomena  that  were  to  be  measured.  The  Jangle- 
HE  series  was  conducted  in  the  fall  of  1951,  and  involved  four  shots: 
three  2,560-pound  charges  and  a 20- ton  event.  Three  of  the  HE  spheres 
were  tangent  to  the  surface  (above  and  below) , and  the  fourth  was  shal- 
lowly buried.  An  additional  objective  of  these  tests  was  to  relate  the 
results  of  the  UET  program,  conducted  at  Dugway,  Utah,  to  the  soil  char- 
acteristics of  the  desert  alluvium  at  NTS. 

Project  Mole  vas  an  HE  program  consisting  of  four  test  series: 
two  in  a sand-gravel  medium,  one  in  dry  clay,  and  one  in  moist  clay. 

All  charges  were  256-pound  TNT  spheres,  with  charge  positions  ranging 
from  6 feet  above  to  6 feet  below  the  ground  surface.  The  purpose  of 
Project  Mole  was  to  investigate  the  relations  between  charge  position 
and  soil  type  and  to  study  the  effects  of  underground  explosions.  The 
program  was  conducted  by  the  Stanford  Research  Institute  for  the  Army 
Corps  of  Engineers  in  19 52-5 

A number  of  studies  were  conducted  at  WES  during  the  period  1957-60 
to  investigate  various  phenomena  associated  with  cratering.  Among 
these  were  investigations  of  the  effects  of  different  types  of  charge 
stemming  (backfill  of  the  charge-emplacement  hole),  the  effects  of  a 
shallow  soil-rock  interface  below  the  explosion  area,  and  the  formation 
of  craters  in  loess  and  clay.  Charge  weights  in  these  tests  ranged  from 
1/6  to  256  pounds. 

2.2  HE  TEST  PROGRAMS:  1962  AND  AFTER 

In  1962,  the  entire  complexion  of  testing  programs  involving  large, 
crater-forming  explosions  changed.  Part  of  the  change  vas  due  to  a 
redirection  of  cratering  research  in  support  of  the  nevly  initiated 
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Plowshare  Program  to  utilize  nuclear  explosives  for  peaceful  purposes 
(discussed  in  Section  2.4),  but  the  increased  intensity  of  testing  ef- 
forts was  in  larger  part  due  to  the  end  of  the  first  moratorium  on 
atmospheric  nuclear  testing  in  1961  and  the  consequent  resumption  of 
nuclear  testing  by  the  U.  S.  early  in  1962.  The  most  important 
cratering  experiments  conducted  during  this  period  are  illustrated  in 
Figure  2.1.  Nuclear  cratering  tests  were  conducted  at  NTS  from  March 
through  July  1962.  After  a six-month  pause,  an  additional  series  of 
experiments,  named  Ferris  Wheel,  was  planned  for  early  1963  to  directly 
compare  the  cratering  efficiencies  of  nuclear  and  TNT  explosives  of 
equal  yield.  The  advent  of  the  second  moratorium  on  nuclear  testing 
came  within  a few  weeks  of  the  execution  of  Ferris  Wheel,  however,  and 
the  program  was  cancelled. 

The  Flat  Top  Series  utilized  two  of  the  experimental  arrays  that 
were  emplaced  in  Frenchmen  Flat  at  NTS  for  the  Ferris  Wheel  Series. 

Event  II  was  a 20-ton,  half-buried,  spherical  TNT  charge  as  originally 
planned  for  Ferris  Wheel,  while  Event  III  involved  the  substitution  of 
an  identical  20-ton  TNT  charge  for  the  originally  planned  nuclear  de- 
vice. A third  event.  Flat  Top  I,  was  added  to  the  series  to  compare  the 
craters  formed  in  playa  silt  by  Events  II  and  III  with  a crater  formed 
by  an  identical  PO-tcn  charge  in  limestone.  The  entire  program  was 
sponsored  and  directed  by  the  Defense  Atomic  Support  Agency,  now  the  De- 
fense Nuclear  Ageacy  (DNA). 

The  Air  Vent  Series  was  sponsored  by  DNA  and  conducted  in  1963-64 
by  the  Sandia  Corporation  to  provide  a tie-in  between  the  Flat  Top 
crater  results  in  playa  and  previous  large  HE  cratering  tests  fired  at 
deeper  depths  of  burst  (DOB)  in  NTS  desert  alluvium.  The  Air  Vent 
Series  was  composed  of  three  phases:  Phase  I was  a single  20-ton  TNT 
sphere  detonated  at  a DOB  of  17  feet,  Phase  II  was  a series  of  tventy 
256-pound  spherical  TNT  charges  fired  at  a vide  range  of  DOB,  and 
Phase  III  was  a scries  of  nine  TNT  charges  of  64,  1,000,  and  6,000 
pounds,  all  fired  as  true  surface  bursts  (center  of  gravity  at  ground 
surface). 

The  Multiple  Threat  Cratering  Experiment  (MICE)  was  sponsored  by 


DIIA  and  conducted  by  the  Boeing  Company  under  the  supervision  of  the 
Air  Force  Weapons  Laboratory  at  the  U.  S.  Army  Yakima  Firing  Range, 

Washington,  in  I965.  The  purposes  were  to  investigate  the  effects  of 
detonating  successive  charges  along  a single  vertical  axis,  and  to 
study  the  influence  of  the  charge  shape  on  the  crater  and  ejecta  from 
near-surface  explosions.  The  t-fTCE  tests  consisted  of  eighteen 
4,000-pound  charges  and  two  16,000-pound  TNT  charges,  with  one  each  of 
the  two  sizes  being  hemispherical  and  the  -emainder  spherical.  Nine  of 

the  spherical  4,000-pound  charges  were  "nail-driving"  shots  (i.e,,  the 
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charge  vas  detonated  in  the  crater  of  a preceding  detonation),  four  were 
true  surface  events,  three  were  surface-tangent  above  the  surface,  and 
one  vas  surface-tangent  below  the  surface.  The  second  16,000-pound 
charge  vas  a true  surface  detonation.  All  events  were  fired  in  a weak 
basalt. 

Additional  nail-driving  experiments  were  conducted  in  small-scale 
experiments  in  1963-64,  prior  to  the  MTCE  tests.  A single  series  of 
64-pound  TNT  charges  vas  fired  in  desert  alluvium  in  a Sandia  Corpora- 
tion test  program,  while  4-  and  21-pound  TNT  charges  were  fired  in  a 
sandy-clayey  silt  by  WES. 

Operation  Sailor  Hat  vas  a series  of  three  500-ton,  hemispherical 
charges  detonated  on  the  surface  of  a basalt  medium  on  Kahoolave  Island, 

Hawaii,  in  1965.  Crater  measurements  vere  only  made  for  the  first  event, 
since  the  latter  tvo  events  vere  decoupled  from  the  basalt  by  placing 
the  charges  on  artificial  fill  materials . Sailor  Hat  vas  conducted  by 
the  Navy,  under  the  sponsorship  of  UNA,  to  test  the  response  of  ships 
(anchored  offshore)  to  airblast  loadings. 

The  Mine  Shaft  Series  vas  conducted  in  1968-69  on  a granite  medium 
near  Cedar  City,  Utah.  These  tests  vere  sponsored  by  DNA- and  conducted 
by  WES  to  develop  data  on  the  effects  of  near-surface  explosions  over  a 
hard  rock  medium.  The  main  events  vere  preceded  by  ten  1,000-pcucd 
calibration  shots  fired  at  different  heights  of  burst  (HOB)  very  near 
the  air-rock  interface.  The  Mine  Shaft  I Series  involved  tvo  100-ton 
spherical  TNT  charges:  Mine  Under  at  an  HOB  of  tvo  charge  radii  (or 

2rch)  and  Mine  Ore  at  an  HOB  of  0.9rfih  . The  Mine  Shaft  II  Series  j 

. 
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consisted  of  the  100- ton  Mineral  Rock  Event,  a duplicate  of  Mine  Ore, 
and  the  Mineral  Lode  Event,  a 16-ton  spherical  charge  of  ammonium 
nitrate  slurry  detonated  at  a DOB  of  100  feet. 

Cratering  experiments  were  conducted  at  the  Canadian  Defence  Re- 
search Board's  Suf field  Experimental  Station  (SES,  later  designated  the 
Defence  Research  Establishment,  Suffield,  or  DRES)  at  Ralston,  Alberta, 
from  1958  to  1970.  Both  the  Watching  Hill  and  Drowning  Ford  test  ranges 
at  DRES  consist  of  glacial  till,  a heterogeneous  mixture  of  clay,  silt, 
gravel,  and  fine  sand  strata.  An  early  te9t  series  involved  the  detona- 
tion of  7**  TNT  charges  ranging  in  weight  from  8 to  10,065  pounds,  in 
both  spherical  and  hemispherical  charge  shapes.  The  spherical  charges 
were  fired  in  both  half-buried  and  surface-tangent  geometries,  while  the 
hemispherical  charges  were  all  resting  on  the  surface.  From  1959 
through  1963,  six  additional  detonations  of  large,  hemispherical  TNT 
charges  occurred  at  DRES,  with  charge  weights  ranging  from  5 to  100  tons. 
In  196b,  a 500-ton  hemispherical  charge  was  detonated  on  the  surface  for 
the  Snowball  Event.  As  with  other  tests  involving  hemispherical  charges, 
airblast  measurements  primarily  dictated  shot  geometry. 

The  Distant  Plain  Series  was  conducted  in  1966-6?  at  DRES  as  a part 
2 

of  the  Quadr apart ite  program  to  develop  Improved  methods  for  simulating 
and  predicting  the  effects  of  nuclear  explosions.  Six  of  the  Distant 
Plain  Events  produced  craters:  Events  1A,  3>  5,  and  6A  were  20-ton 
spherical  TNT  charges,  with  the  first  fired  at  an  HOB  of  29  feet,  the 
second  and  third  half -buried,  and  the  last  placed  tangent  to  the  sur- 
face; Event  4 was  a 50-ton  hemispherical  charge  detonated  on  the  ground 
surface  in  a forest  near  Hinton,  Alberta;  and  Event  6 was  a 100-ton 
spherical  charge  in  a surface-tangent  geometry. 

Two  recent  tests  were  made  at  DRES,  both  500-ton  spherical  charges 
fired  surface-tangent.  The  Prairie  Flat  Event  occurred  in  1968  and  the 
Dial  Pack  Event  in  1970. 


O 

A cooperative  effort  for  blast-effects  research  between  the  U.  S. , 
Great  Britain,  Canada,  and  Australia.  Prior  to  the  inclusion  of 
Australia,  it  was  known  as  the  Tripartite  program. 
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2.3  NUCLEAR  TESTS 

2.3.1  Pacific  Test  Programs.  The  testing  of  nuclear  devices  for 
weapons  development  began  at  the  Pacific  Test  Range  almost  immediately 
after  World  War  11.  Of  the  many  testa  conducted,  only  a few  produced 
measurable  craters  (i.e.,  the  major  portion  of  the  crater  being  con- 
tained in  a land  area),  and  relatively  few  of  those  craters  were  ade- 
quately measured.  Ten  of  the  Pacific  tests  produced  craters  whose 
measurements  were  recorded.  All  were  fired  at  either  the  Bikini  or 
Eniwetok  Atolls,  vhere  the  soil  medium  is  described  as  coral  sand  under- 
lain with  intermittent  beds  of  hard  coral  and  cemented  rubble. 

The  Mike  Event  of  Operation  Ivy  in  1952  was  the  first  te3t  of  a 
thermonuclear  weapon.  The  10.U-Mt  device  was  fired  at  a height  of 
35  feet  above  the  island  surface.  Five  more  high-yield  devices,  ranging 
from  1.3  to  15  Mt,  were  detonated:  Bravo,  Zuni,  and  Tewa  at  Bikini  and 

Koa  and  Oak  at  Eniwetok.  Tewa  and  Oak  were  both  detonated  on  barge3 
floating  in  shallow  water.  The  Koon  Event  was  a 110-kt  detonation  at 
Bikini,  while  Lacrosse,  Seminole,  and  Cactus  were  U0-,  1U-,  and  18-kt 
yields,  respectively,  at  Eniwetok.  Seminole  was  detonated  inside  a tank 
of  water. 

2.3.2  NTS  Te3t  Programs.  The  first  nuclear  cratering  experiments 
conducted  at  the  NTS  were  the  Jangle-S  and  Jangle-U  Events  of  Operation 
Plumbob  in  1957.  Both  were  1.2-kt  yields,  fired  in  desert  alluvium — 
Jangle-S  at  an  HOB  of  3.5  feet  and  Jangle- J at  a DOB  of  17  feet.  The 
Teapot  Ess  Event  followed  in  1955,  and  was  detonated  still  deeper  in 
alluvium  at  67  feet.  In  1958,  the  Neptune  Event  was  fired  inside  a 
chamber  100  feet  beneath  a sloping  mountainside,  in  a volcanic  tuff 
medium.  The  Neptune  device  had  a yield  of  0.115  kt. 

No  nuclear  cratering  tests  were  conducted  during  the  first  mora- 
torium on  nuclear  testing  from  1959  through  1961.  After  the  resumption 
of  testing  by  the  USSR  in  1961,  the  0.^2-kt  Danny  Boy  Event  was  fired 
in  March  19&2,  at  a DOB  of  110  feet  iu  a basalt  mesa  at  NTS.  Danny  Boy 
was  soon  followed  by  Small  Boy,  a low  airburstJ  over  playa,  and  the 

3 

Classified  yields  are  omitted. 


Little  Feller  I and  Little  Feller  II  Event3,  all  in  desert  alluvium. 
Johnie  Boy  was  a near-surface  detonation,  while  the  Little  Feller  Events 
were  both  small  weapons  fired  at  low  HDB's.  Except  for  the  Plowshare 
Program,  no  further  nuclear  cratering  tests  were  scheduled  until  the 
Ferris  Wheel  Events  of  February  1963,  which  were  cancelled  due  to  the 
commencement  of  the  second  moratorium  on  nuclear  testing. 

2.4  PLOWSHARE  EXPERIMENTS 

The  Plowshare  Program  was  inaugurated  by  the  U.  S.  Atomic  Energy 
Commission  (AEC)  in  1957  to  develop  peaceful  applications  of  nuclear 
explosive  energy.  The  most  immediately  obvious  application  was,  of 
course,  earthmoving  by  explosive  cratering.  The  first  actual  planned 
use  of  nuclear  cratering  for  constructive  purposes  was  the  excavation 
of  a harbor  on  the  coast  of  Alaska.  This  experiment,  nicknamed  Project 
Chariot,  was  to  proof-teat  the  feasibility  of  nuclear  excavation  in 
the  early  1960's.  Unfortunately,  the  project  was  soon  postponed,  then 
eventually  cancelled.  Among  the  planned  or  proposed  projects  for 
nuclear  excavation  were  (l)  the  excavation  of  a railroad  pass  through 
a mountain  ridge  in  California,  (2)  the  construction  of  a portion  of 
a canal  connecting  the  Tennessee  and  Tombigbee  Rivers,  (3)  the  olasting 
cf  a harbor  at  Cape  Keraudren,  Australia,  and  (4)  a sea-level  canal 
through  or  near  the  Isthmus  of  Panama.  As  of  this  writing,  the  use  of 
nuclear  explosions  as  & means  of  producing  excavations  has  been  indef- 
initely postponed  as  either  unjustifiable,  uneconomical,  or  unpopular 
from  an  environmental  and/or  political  viewpoint. 

The  most  ambitious  concept  for  nuclear  excavation  was  the  conntruc- 
tion  of  a new,  sea-level  "Panama"  Canal.  This  project  afforded  the  pri- 
mary stimulus  for  extensive  research  in  the  field  of  cratering  technology 
supported  by  the  AEC  for  the  past  decade  under  the  Plowshare  Program. 

Much  of  the  testing  was  done  in  desert  alluvium  at  NTS  because  of  its 
availability  and  the  economic  advantages  of  support  at  NTS,  but  many  of 
media  later  selected  for  cratering  experiments  were  chosen  due  to  their 
similarity  to  the  soil  and  rock  indigenous  to  the  Panamanian  Isthmus. 

The  first  significant  cratering  experiment  under  the  Plowshare 


Program  was  Project  Stagecoach  which  included  three  separate  detonations 
of  20-ton  spherical  TNT  charges  at  DOB's  of  17,  34,  and  30  feet.  The 
tests  were  conducted  in  desert  alluvium  at  NTS  in  March  1950.  One  of 
the  chief  objectives  of  Stagecoach  wa3  to  determine  if  cube-root  scaling 
was  actually  valid  for  HE  crater  formation  in  desert,  alluvium.  Sandia 
Corporation  was  primarily  responsible  for  the  conduct  of  this  experiment 
and  the  next  two  which  followed  it. 

Project  Buckboard  was  a similar  test  aeries  fired  in  basalt  at  the 
NTS  to  provide  basic  information  on  the  formation  of  craters  in  a hard 
rock  environment.  This  series,  conducted  in  I960,  consisted  of  ten 
1,000-pound  cast  TNT  charges  fired  at  DOB's  ranging  from  5 to  25  feet  in 
5-foot  intervals,  followed  by  three  40,000-pound  spherical  TNT  charges 
fired  at  DOB’s  of  26,  43,  and  60  feet. 

Project  Scooter  was  a 500-ton  spherical  TNT  charge  fired  at  a DOB 
of  125  feet  in  desert  alluvium  at  NTS  in  October  I960.  This  test,  the 
final  in  the  Sandia  series,  was  designed  to  extend  knowledge  of  the 
mechanics  of  crater  formation  into  yields  near  the  kiloton  regime. 

The  Sedan  Event  was  the  largest  explosive  cratering  experiment  ever 
conducted  by  the  U.  S.  The  100-kt  device  was  detonated  in  1962  at  a DOB 
of  635  feet  in  tho  desert  alluvium  at  MTS.  The  main  purpose  of  the  ex- 
periment was  to  extend  empirically  based  cratering  theory  to  large  yields 
representative  of  those  that  would  be  employed  in  an  actual  large-scale 
construction  project  under  the  Plowshare  Program.  Equally  important  was 
the  evaluation  of  the  physical  hazards,  particularly  radiation,  that 
could  be  expected  frco  such  employment.  The  Lawrence  Radiation  Labora- 
tory, now  the  Lawrence  Livermore  Laboratory  (LLL)  at  Livermore,  Cali- 
fornia, had  primary  responsibility  for  the  Sedan  experiment. 

In  late  1962,  the  Nuclear  Cratering  Group,  now  the  WES  Explosive 
Excavation  Research  Laboratory  (EERL),  was  created  as  an  organization 
under  the  U.  S.  Army  Corps  of  Engineers  to  have  direct  responsibility 
for  supervising  experimental  research  on  the  use  of  nuclear  explosions 
for  civil  construction  purposes.  The  first  test  program  under  the  su- 
pervision of  EERL  was  Pre-Buggy  I,  conducted  in  desert  alluvium  at  NTS 
in  the  winter  of  1962-63  to  develop  design  criteria  for  a future 
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large-scale  nuclear  row-charge  cratering  experiment,  Pre-Buggy  I was,  in 
turn,  parcially  designed  from  data  acquired  in  an  earlier,  smaller  scale 
row  experiment  called  Froject  Rowboat  (by  LLL).  Six  single-charge  tests 
were  fired  to  determine  optimum  DOB,  followed  by  four  5-charge  row 
events  to  determine  the  most  desirable  charge  spacing  for  the  creation 
of  row  craters  or  channels.  Each  charge  wa3  a 1,000-pound  sphere  of  a 
liquid  explosive,  nitromethane  (NM).  The  Pre-Buggy  II  Series  was  fired 
several  months  later,  and  consisted  of  six  row  shots  in  which  charge 
spacings,  DOB,  and  stemming  were  further  varied. 

The  Pre-Schooner  Series  wa3  executed  by  EERL  in  Februury  1964  in 
basalt  on  the  Buckboard  Mesa  at  NTS.  Four  shots  of  20  tons  of  NM  each 
were  fired  at  DOB's  ranging  from  42  to  66  feet.  These  were  soon  fol- 
lowed by  the  Dugout  Event , a row-charge  crater  shot  in  basalt  consisting 
of  five  20-ton  NM  charges.  In  December  of  196k,  a low-yield  (85-ton) 
nuclear  cratering  event  named  Sulky  was  detonated  in  the  basalt  at 
Buckboard  Mesa.  The  DOB  of  90  feet  was  slightly  greater  than  optimum 
in  an  attempt  to  contain  a greater  percentage  of  radioactivity.  The 
Pre-Schooner  II  Event  was  an  85 -ton  charge  of  NM  fired  in  rhyolite  in 
southwestern  Idaho  in  September  1965. 

The  Palanquin  Event  was  a 4.3-kt  nuclear  cratering  experiment  in  a 
rhyolitic  rock  at  NTS,  fired  at  a DOB  df  280  feet,  in  April  1965.  A 
second  nuclear  cratering  experiment  in  rhyolite,  the  2.3-kt  Cabriolet 
Event,  was  fired  in  January  1968.  The  largest  cratering  shot  in  rock 
was  the  35-kt  Schooner  Event,  fired  in  tuff  at  NTS  in  December  1968.  In 
1969  the  Buggy  Event,  a nuclear'  row  charge  consisting  of  five  1.1-kt  de- 
vices, was  detonated  in  basalt  at  NTS. 

An  extensive  program  of  testing  known  as  Pre-Gondola  was  also  con- 
ducted by  EERL  in  a wet  clay-shale  medium  at  Ft.  Peck  Reservoir,  Montana, 
from  October  1966  to  October  1968.  Charge  weights  ranged  fraa  very  small 
(less  than  10  pounds)  to  40  tons,  in  both  single-charge  crater  tests  and 
various  row-charge  tests.  Most  of  these  experiments  employed  NM  as  the 
explosive,  although  sane  tests  used  an  ammonium  nitrate  slurry. 
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Explosive  equivalences  in  terms  of  cratering  efficiency  are  considered 
in  Chapter  3. 
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Concurrent  with  their  larger  scale  experiments,  EERL  carried  out 
an  extended  series  of  3ingle-  and  row-charge  studies  using  1-pound 
charges  in  a closely  controlled  sand  medium,  under  the  name  Project  Zulu. 

An  additional  program  of  multiple-charge  cratering  research,  using 
relatively  small  HE  charges,  has  been  conducted  at  Sandia  for  the  last 
decade  under  sponsorship  of  the  Plowshare  Program.  Most  of  these 
studies  have  been  concerned  with  the  cratering  effects  of  multiple 
charges  detonated  simultaneously  or  in  a particular  sequence,  such  as 
adjacent  rows  of  charges  fired  in  sequence,  square  arrays  of  charges 
fired  simultaneously,  or  vertical  arrays  of  charges.  Virtually  all  of 
these  tests  were  conducted  in  desert  alluvium  near  Albuquerque,  New 
Mexico. 

2.5  USSR  CRATERING  EXPERIMENTS 

In  recent  years,  an  ambitious  program  of  civil  cratering  applica- 
tions has  been  undertaken  in  Russia,  generally  paralleling  the  U.  S. 
Plowshare  Program.  At  this  writing,  11  nuclear  projects  have  been  car- 
ried out,  as  described  in  Reference  3.  Since  reported  results  are  in- 
complete, these  tests  are,  with  one  exception,  not  included  in  this 
report.  Reports  of  these  projects  mAke  interesting  reading,  however, 
and  they  are  listed  in  Appendix  C. 

2.6  SUMMARY 

The  foregoing  synopsis,  which  lists  only  major  recent  cratering 
experimentation  (illustrated  in  Figure  2.1),  indicates  the  accumulation 
of  a wealth  of  crater  data,  much  of  it  in  the  past  decade.  To  the 
casual  reader,  it  might  appear  that  sufficient  testing  has  been  accom- 
plished to  answer  virtually  all  questions  pertaining  to  this  subject; 
unfortunately,  this  is  not  the  case.  Cratering  phenomena  are  so  com- 
plex, snd  cratering  applications  are  so  varied,  that  much  yet  remains 
to  be  done.  In  the  chapters  that  follow,  an  attempt  will  be  made  to 
define  basic  crater  parameters,  drawing  on  all  applicable  data  available 
to  the  authors,  and  methods  of  applying  this  information  will  be  dis- 
cussed. In  the  process,  areas  of  insufficient  data  will  be  pointed  out. 
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CALENDAR  YEARS 


Illustrated  history  of  major  blast-  and  shock-effects  tests 


CHAPTER  3 


CRATERING  DATA 


3.1  DESCRIPTION  OF  DATA 

The  site  and  shape  of  an  explosively  generated  crater  are  dependent 
upon  many  different  factors.  The  areas  of  influence  are  divided  into 
two  main  categories:  environmental  parameters  and  explosive  parameters. 

The  environmental  parameters  include  the  density,  composition,  and 
strength  of  the  medium  in  which  the  crater  is  formed,  the  material's 
moisture  content,  the  Jointing  or  layering  of  the  medium,  and,  for  the 
sake  of  completeness,  the  atmospheric  conditions  at  the  time  of  detona- 
tion. Explosive  parameters  include  the  quantity  and  type  of  explosive 
used,  the  geometry  of  the  charge,  the  method  of  charge  emplacement , 
and  the  position  of  the  charge  relative  to  the  air-medium  interface, 

(in  every  case,  a complete  detonation  of  the  charge  is  assumed).  The 
actions  and  interactions  of  these  independent  variables  on  crater  forma- 
tion are  complex,  and  only  the  major  influencing  factors  are  examined 
in  this  report.  These  major  factors  form  the  basis  upon  which  the  data 
tables  (Appendix  A)  are  subdivided.  Within  the  appropriate  subdivisions , 
data  on  crater  sixes  and  shapes  are  included. 

The  tables  are  as  complete  a3  is  possible  in  a mass  compilation  of 
this  sort.  If  all  variables  have  tabulated  values  for  any  particular 
shot,  the  crater  can  be  essentially  reconstructed.  For  example,  by 
using  the  available  apparent  crater  depth,  radius,  lip  height,  and  angle 
of  slope  in  connection  with  a crater  shape  code,  a typical  crater  pro- 
file can  be  drawn. 

3.2  EXPLOSIVES  CHARACTERISTICS 

Approximately  10  different  kinds  of  explosive  have  been  widely  used 
for  different  test  programs,  depending  upon  the  requirements  of  the  test 
and/or  the  convenient  availability  of  various  explosive  types.  For  this 
report,  all  explosives  are  converted  to  an  equi valent  TNT  yield*  each 
equivalent  TNT  weight  is  assumed  to  be  a sphere  of  cast  explosive. 
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The  explosives  that  have  been  most  commonly  used  for  cratering 
research  are  TNT,  NM,  ammonium  nitrate/fuel  oil  mixtures  (ANFO),  and 
ammonium  nitrate  (AN)  slurries.  TNT  has  been  selected  as  the  basic  ex- 
plosive due  to  its  widespread  use  in  the  past,  a result  of  the  fact  that 
it  is  readily  available  (to  military  users);  it  can  be  cast  in  any  size 
and  shape,  and  it  is  a relatively  insensitive  explosive  and  thus  safe 
to  handle.  NM  is  a liquid;  actually,  it  is  classed  as  a solvent  rather 
than  an  explosive.  It  has  been  widely  used  in  Plowshare  cratering  ex- 
periments due  to  the  fact  that  a large  volume  of  NM  (up  to  100  tons  has 
been  used)  can  be  poured  into  an  underground  emplacement  chamber  through 
small-diameter  fill  pipes  from  the  surface,  thus  greatly  simplifying  the 
problem  of  charge  emplacement  for  large,  deeply  buried  shots.  ANFO  has 
been  used  for  many  years  for  cratering.  It  is  classed  as  a blasting 
agent,  and  can  be  mixed  on-site  prior  to  placement  in  its  cavity  or  con- 
tainer. AN  slurry  has  been  used  increasingly  in  recent  years,  since  it 
is  also  classed  as  a blasting  agent  and  can  be  pumped  into  deeply  buried 
charge  emplacement  chambers;  its  higher  viscosity  in  most  instances 
eliminates  the  need  for  extensive  lining  of  the  chamber  to  prevent 
leakage . 

Unfortunately,  very  few  data  exist  that  permit  performance  compari- 
sons of  the  crater-forming  abilities  of  different  types  of  chemical  ex- 
plosives as  compared  to  the  base  explosive  TNT.  In  the  absence  of  more 
definitive  research  on  cratering  performance,  a list  of  common  explo- 
sives and  their  TNT  equivalents  has  been  compiled  from  References  1 and 
10  and  is  presented  as  Table  3*1.  Table  3.1  does  not  take  into  account 
the  effects  of  different  soils  or  rocks  on  cratering. 

3.3  DATA  PRESENTATION 

All  crater  data  available  to  the  authors  and  falling  within  the 
scope  of  this  report  are  presented  in  Appendix  A,  Tables  A.l  through 
A. lb.  The  data  tables  are  subdivided  into  three  major  parts.  First, 
they  are  separated  according  to  type  of  medium,  since  crater  size  is 
strongly  dependent  on  strength  of  the  material  cratered.  Second,  the 
shots  in  each  medium  are  grouped  according  to  water  content  because  of 
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the  variation  in  medium  strength  and  shock-propagation  characteristics 
with  change  in  water  content.  For  example,  all  crater  data  for  sand 
are  further  subdivided  into  wet,  dry-to-moist , and  dry  sands.  The  third 
data  group  consists  of  10  categories  of  scaled  depths  (heights)  of  burst, 
since  differences  in  charge  positioning  affect  energy  partitioning  into 
the  medium.  For  specifying  the  charge  position,  cube-root  scaling  has 
been  used.  Thus,  depth  or  height  of  the  charge  position  is  divided  by 
the  cube  root  of  charge  weight,  permitting  the  grouping  of  explosion 
events  of  different  yields.  Selection  of  this  cube-root  scaling  expo- 
nent is  further  discussed  in  Chapter  4.  Note  that  the  categorical  divi- 
sions are  uneven;  they  have  been  selected  to  identify  HOB/dOB  at  which 
critical  changes  in  crater  dimensions  and/or  crater  scaling  "laws"  can 
be  expected. 

A special  data  group  which  stands  apart  from  the  regular,  spherical- 
charge  cratering  tests,  but  which  is  properly  included  in  a cratering 
guide  such  as  this,  is  that  of  the  hemispherical-charge  data  contained 
in  Table  A.14.  All  available  cratering  data  for  this  charge  geometry, 
which  is  frequently  encountered  in  connection  with  airblast  measurements , 
are  included  in  this  table. 

The  data  in  Tables  A.l  through  A. 13  are  graphically  presented  in 
Appendix  B,  Figures  B.l  through  B-92;  Figure  B.93  is  a graphical  repre- 
sentation of  the  hemispherical-charge  data  from  Table  A.l4.  Data  are 
fitted  by  the  method  of  least  squares  to  provide  empirical  equations. 

Each  equation  is  based  on  a minisum  of  two  data  points  representing  a 
minimum  range  in  charge  yields  *2/3  of  a logarithmic  cycle.  In  addition, 
no  curve  with  a slope  1.5  < ^ < 0 was  considered.  Each  curve  was 
drawn  from  1 pound  on  the  abscissa  to  a point  slightly  beyond  the  largest 
yield  for  which  data  were  available. 

The  computer  program  used  to  separate,  analyte,  and  plot  the  data 
is  shewn  in  Appendix  D.  Explanations  necessary  to  the  use  of  the 
crater-data  tables  precede  the  tables.  Bote  especially  -that  only  HOB 
is  listed;  a negative  number  denotes  a buried  charge  geometry  (DOB). 
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3.U  SUBSURFACE  DEFORMATION 


J 


As  shown  in  Figure  1.2,  extensive  permanent  subsurface  deformation 
accompanies  a large  crater.  The  outer  limit  of  this  deformation,  too 
indistinct  for  precise  measurement,  marks  the  boundary  at  which  the 
transient  -tresses  (compressive,  tensile,  or  shear)  imposed  by  the  ex- 
plosion shock  wave  exceed  the  respective  medium  strengths.  Beyond  this 
boundary,  the  medium  responds  elastically. 

Subsurface  deformation  may  be  divided  into  two  general  zones: 
rupture  and  plastic  deformation.  The  rupture  zone  is  so  named  because 
of  the  extensive  crushing  and  cracking  which  it  contains,  while  the 
plastic  zone  exhibits  a smoother  flovage  of  the  medium  (in  soil),  often 
unnoticeable  unless  special  methods , such  as  preemplaced  colored  sand 
columns,  are  used  to  detect  such  movement.  In  the  upper  portion  of  the 
plastic  zone,  a shearing  action  often  occurs  by  slippage  along  hori- 
zontal planes  of  weakness.  In  rock,  the  plastic  zone  is  nonexistent 
or  insignificant,  but  permanent  displacement  occurs  by  the  closing  or 
opening  of  Joints.  In  rock  as  well  as  many  soils,  boundaries  of  the 
individual  zones  of  deformation  are  irregular  and  indistinct. 

Relatively  few  observations  are  available  for  definition  of  sub- 
surface deformation  zones,  especially  for  large  explosions.  Summarized 
in  the  following  tabulation,  in  terns  of  true  crater  radius  r^  and 
depth  d.  or  cavity  radius  r , are  "rule-of-thumb"  limits  of  these 
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zones  in  soil  and  rock: 


0 


Type  of 

Shot  Geometry 

r flUlui  " t 

Near-Surface 

Deeply  Buried 

Horizontal  Vertical 
Limit  Limit 

Horizontal 

Limit 

Vertical 

Limit 

For  Soil 

Rupture 

2*t 

2rt 

1.54,  * 

Plastic 

6rt  3.04, 

(Continued) 

\ 

__b 

Where 


d m DOB  ♦ r . 
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Insufficient  data. 
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Type  of 
Failure 

Shot  Geometry 

Near-Surface 

Deeply  Buried 

Horizontal 

Limit 

Vertical 

Limit 

Horizontal 

Limit 

Vertical 

Limit 

For 

Rock 

Rupture 

2rt 

3.0dt 

2rt 

2.0dta 

Displacement 

4rt 

4.0dt 

b 

b 

* Where  <L  ■ DOB  + r . 

v t C 

Insufficient  data. 


Dimensions  of  cavity  radius  r are  discussed  in  Section  6.U. 

c 

Plastic  deformations  or  displacements  are  intended  to  include  only 
those  which  are  definitely  measurable.  In  general,  movement  of  the 
medium  close  to  the  charge  is  radially  outward,  or  away  from  the  explo- 
sion. There  is  often  a vertical  component  to  such  movement,  usually 
upward  for  surface  and  buried  charges  (upthrust)  and  downward  for  above- 
surface geometries.  Farther  from  the  charge,  permanent  displacement 
toward  GZ  has  frequently  been  observed.  Some  surface  rotation  (in  the 
horizontal  plane)  is  usually  observed  on  survey  monuments  that  were  per- 
manently displaced  by  an  explosion,  but  no  pattern  of  rotation  is  dis- 
cernible. Relaxation  of  the  deformed  medium  has  also  been  observed  to 
occur  over  a period  of  several  days,  while  a condition  of  equilibrium  is 
established  among  pore  and/or  Joint  pressures  and  the  crater  void  itself. 
Relaxation  is  manifested  by  a small  but  measurable  movement  of  the  com- 
pacted medium  tcward  GZ. 
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TABLE  3.1  COMPARISON  OF  EXPLOSIVE  CRATERING  EFFICIENCY 
WITH  THAT  OF  TNT  ( REFERENCES  1 AND  10) 


) 


To  determine  relative  cratering  efficiency  (TOT) , mul- 
tiply weight  of  explosive  charge  by  conversion  factor. 


Explosive 

Conversion  Factor 

TNT* 

1.00 

Amatol 

0.91* 

Dynamite  (koj) 

0.68 

Pentolite 

1.23 

C-k,  C-3 

1.3k 

Ammonium  Nitrate 

1.00 

Nitromethane 

1.10 

* TNT  explosive  energy  **  10^  calorie s/ton  heat  of  deto- 
nation * 1012  calories/kt. 
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CHAPTER  4 

ANALYSIS  OF  DATA  TRENDS 

4.1  OUTER  SIZE  AND  SHAPE 

The  logarithmic  graphs  of  Appendix  B are  summarized  in  Figures  4.1 
through  It. 19.  Figures  4.1  through  U .9  are  constructed  to  show  crater  di- 
mensions at  the  1-ton  yield  (weight)  level  in  terms  of  TNT  equivalence 
for  those  media  for  which  sufficient  data  exist  and  for  the  range  of 
HOB/DOB  available.  For  this  yield,  crater  radius,  depth,  and  lip  height 
may  be  read  directly  from  the  curves;  for  other  yields,  the  graphs  may 
be  entered  with  scaled  HOB  or  DOB,  and  dimensions  cay  be  scaled  to  the 
appropriate  charge  weight  using  the  nearest  value  of  the  scaling  expo- 
nent shown  or  an  interpolated  value.  Section  4.3  discusses  scaling  In 
more  detail.  Figures  4.10,  4.11,  and  4.12  each  show  the  variation  of 
a single  crater  dimension  for  the  same  charge  weight  with  different 
media.  This  makes  possible  an  estimate  of  crater  dimensions  for  those 
media  for  which  insufficient  data  exist  to  support  a separate  curve. 

' Figures  4.13  through  4.19  show  crater  dimensions  for  the  1-kt  yield 
level.  It  will  be  noted  that  fewer  data  are  available  for  this  yield, 
a problem  which  worsens  rapidly  as  yield  increases.  No  attempt  is  made 
in  this  study  to  illustrate  crater  dimensions  above  the  1-kt  level;  the 
general  trends  in  yield  scaling  exponents  for  crater  radius  and  depth 
from  small  to  very  large  yields  are  illustrated,  however,  in  Figure  4.20. 
Scaled  values  for  yields  of  interest  may  be  approximated  by  the  use  of 
Figure  4.21. 

An  Illustration  of  the  geueral  shapes  of  craters  for  charges  of 
all  sizes  and  for  geometries  ranging  from  low  airbursts  to  containment 
depths,  at  which  subsidence  craters  may  be  formed,  is  shown  in  Fig- 
ure 4.22.  Details  of  crater  nomenclature  were  shown  in  Figure  1.2. 
Deeply  buried  explosions  are  discussed  in  more  detail  in  Section  6.4. 

4.2  CRATERING  MECHANISMS 

An  understanding  of  the  mechanics  of  crater  formation  aids  in 


C 


37 


predioting  (or  explaining)  the  size  and  shape  of  the  crater  void  and  the 
volume  of  material  ejected  from  it.  In  addition  to  the  three  primary 
cratering  variables  (charge  yield,  shot  geometry,  and  characteristics 
of  the  cratered  medium) , there  are  three  basic  mechanisms  which  govern 
the  formation  of  HE  craters:  material  ejection,  compaction,  and  plasr 
tic  flowage.  For  practical  purposes,  the  same  may  be  said  of  nuclear 
craters,  since  the  volume  of  material  involved  in  the  fourth 
mechanism — vaporization— is  generally  regarded  as  being  insignificant. 

One  method  of  defining  the  importance  of  different  cratering  mech- 
anisms is  through  their  contributions  to  crater  volume.  A number  of 
volumetric  parameters  for  craters  have  been  isolated  and  defined.  For 
craters  in  soil,  the  following  relations  are  valid: 


vt  * Ydis  + vc  + vf 

(4.1) 

h.  * , 

vdis  " T {ve  + Vfb} 
0 

(4.2) 

Y_  * V.  - V 

fb  t a 

(4.3) 

v,  ■ Kv  + v 

(4.4) 

feu 

V » V, 

u f 

(4.5) 

Where:  v^  » volume  of  true  crater 

v^g  ■ volume  (preshot)  of  material  dissociated  by  the  explosion 

v ■ volume  of  crater  due  to  compression  (compaction) 
c 

vf  ■ volume  of  crater  due  to  plastic  flowage  of  the  medium 

7 >7^  ■ preshot  (in  situ)  and  postshot  unit  weights  of  cratered 
material,  respectively 

v » volume  of  ejecta 

C 

v_  ■ volume  of  fallback 
fb 

v » volume  of  apparent  crater 

v^  ■ total  volume  of  crater  lip 

K » a constant  representing  the  fraction  of  ejecta  volume 
contributing  to  the  formation  of  the  crater  lip 
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vu  ■ volume  of  upthrust  region 

Equations  4.1  through  4.5  permit  assessment  of  the  contributions  of 
various  cratering  mechanisms  to  total  crater  volume.  Similar  relations 
have  been  developed  on  the  basis  of  mass  (Reference  4).  These  same 
equations  (4.1  through  4.5)  are  equally  applicable  in  rock,  except  that 
the  plastic  flowage  term  for  this  medium  is  negligible;  therefore,  Equa- 
tion 4.5  does  not  apply.  Upthrusting  of  the  lip  region  does  occur, 
however,  by  upward  displacement  of  rock  strata  (buried  shots)  and  ex- 
pension  of  joints  by  rebound  action  in  near-surface  shots.  A similar 
action  has  been  noted  beneath  the  rock  craters  for  near-surface  geome- 
tries. When  this  occurs, 


vt  " vdis  - vexp 


(4.6) 


where  v ■ net  volume  of  Joint  expansion, 
exp 

Reference  5 examines  in  detail  the  dimensional  analyses  which 
serve  as  background  to  scaling  relations.  Selecting  from  the  listed 
properties  of  medium  variables  those  which  are  probable  primary  con- 
tributors to  the  three  basic  mechanisms,  the  following  tabulation  can 
be  made: 


Cratering 

Contributing 

Mechanism 

Properties 
of  the  Medium 

Material  Ejection 

c 

Compaction 

vr 

Plastic  Flow 

a,  v 

Where:  o * compression,  shear,  and  tensile  strengths  or  clastic  prop- 

erties ML"*T“‘*  in  units  of  masa-length-time 
v_  ■ void  ratio,  dimensionless 

r _i  _i 

v ■ dynamic  viscosity,  ML  T 


4.3  SCALING  CONSIDERATIONS 


4.3.1  Scaling  aa  a Prediction  Tool.  A primary  purpose  of  the 
graphs  in  Appendix  B is  the  development  of  exponents  to  permit  pre- 
diction of  crater  size.  Generally  referred  to  as  "scaling,”  this  is  a 
mathematical  exercise  of  the  form 

r,  d,  v * kW*  (1».T) 

Where:  r ■ crater  radius  (apparent  or  true) 

d « crater  depth 
v ■ crater  volume 

k * a constant  numerically  equal  to  the  graphical  intercept  at 
V ■ 1 

W ■ charge  weight  (TNT  equivalent) 
n ■ the  scaling  exponent  (slope  of  logarithmic  graph) 

A related  procedure  can  he  employed  when  crater  dimensions  from  a 
similar  experiment  are  available  and  when  u is  known  or  can  be 
approximated.  Thus, 


where  the  subscripts  1 and  2 refer  to  two  different  cratering  events, 
the  results  of  which  are  known  for  one  or  the  other.  The  foregoing 
equations  have  been  used  extensively. 

An  alternative  method,  now  under  development,  is  computer  simula- 
tion of  cratering  processes.  A thorough  discussion  of  this  technique 
is  beyond  the  scope  of  this  report;  briefly,  it  utilizes  a numerical 
computer  routine  such  as  a finite  element  program  to  simulate  the  effects 
of  a detonation  upon  a medium  of  known  physical  properties.  The  results 
can  be  graphed  to  illustrate  crater  formation  at  any  stage  (as  a function 
of  time).  Obviously,  a detailed  knowledge  of  both  the  explosive  and  the 
medium  characteristics  is  prerequisite.  While  this  method  has  contrib- 
uted greatly  to  the  science  of  cratering,  it  will  probably  never  replace 
scaling  as  a rapid  means  of  crater  prediction. 
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4.3.2  Linear  Scaling  Relations.  Reference  5 develops  four  scaling 
rules  governing  crater  radius  and  depth  dimensions.  Each  rule  is  pred- 
icated upon  certain  assumptions,  some  controversial.  Three  of  these 
rules  result  in  the  generally  well-known  cube-root  scaling  (W*^);  the 
fourth  rule,  based  on  energy-gravity  scaling  (charge  energy  and  gravi- 
tational acceleration  g are  considered  when  forming  the  dimensional 
terms,  but  with  g and  the  medium  density  p held  constant),  results  in 
fourth-root  scaling  of  linear  dimensions,  in  all  cases,  experimental 
similarity — to  include  scaling  of  material  properties — is  necessary 
for  unqualified  applicatic--.  of  the  rules.  In  practice,  similarity  re- 
quirements are  seldom  (if  ever)  met,  and  violation  of  these  requirements 
generally  results  in  larger  crater  dimensions  for  Increased  charge  sizes 
than  would  be  predicted  by  formal  scaling  rules.  Thus,  observations  of 
crater  data  frequently  show  n > l/3  , especially  for  lew  yields.  By 
the  same  token,  where  experiment  shews  n < l/3  , the  influence  of 
energy-gravity  scaling  considerations  may  be  suspected.  Since  empirical 
exponents  for  linear  crater  dimensions  show  a general  decrease  with  in- 
creasing charge  yield  (Figure  4.20),  it  seems  likely  that  energy  release 
rather  than  charge  mass  or  weight  becomes  increasingly  important  in  the 
larger  yields.  It  can  be  seen  from  the  foregoing  discussion  that  pre- 
dictions which  span  a wide  difference  in  charge  yields  may  be  subject  to 
the  conflicting  influences  of  similarity  violations  and  energy-gravity 
considerations . 

As  explosion  yield  increases,  and  one  moves  out  of  the  HE  domain 
and  into  the  HE  <u,rain,  basic  differences  in  scaling  behavior  between 
the  two  must  be  considered.  HE  occupies  some  finite  volume  and  gener- 
ates its  own  explosion  gases,  while  HE  is  essentially  a point  source  of 
thermal  energy  which  vaporizes  the  adjacent  material,  thereby  generating 
gases.  Thus,  differences  in  energy  partitioning  and  coupling  into  the 
cratered  medium  almost  surely  cause  differences  in  the  means  by  which 
the  craters  are  formed  and  in  the  scaling  factors  applied  to  them.  Ex- 
perience shows  that  HE  craters  are  larger  than  those  for  comparable  HE 
yields. 

It  seems  probable  that  no  single  scaling  exponent  will  ever 


suffice  to  precisely  predict  any  crater  dimension  except  under  closely 
controlled  conditions  and  for  a limited  range  of  charge  yields.  This  is 
probably  true  in  homogeneous  media,  not  to  mention  real-world  soils, 
where  layering  and  inhomogeneities  abound.  It  is  for  this  reason  that 
empirical  approaches  to  crater  scaling  are  so  generally  used.  Analyses 
presented  in  this  report  show  that  scaling  exponents  for  linear  crater 
dimensions  fall  roughly  between  0.20  and  0.40  for  spherical  charges; 
however,  higher  values  have  been  observed  for  scaling  crater  depths  from 
hemispherical  charges  in  plastic  soil. 

Presumably,  scaling  of  crater  ejecta  field  dimensions  could  be  ac- 
complished in  a manner  similar  to  that  discussed  here.  As  with  crater 
formation,  attempts  are  under  way  to  describe  ejecta  deposition  by  com- 
puter simulation,  as  well  as  by  analysis  of  experimental  data.  Thus  far, 
rather  tenuous  scaling  exponents  in  the  range  of  |o.3  > n > 0.1T  have 
been  developed  from  the  latter,  with  the  larger  exponent  applying  to 
surface  geometries.  Empirical  scaling  exponents  for  ejecta  ranges  are 
further  discussed  in  Chapter  5. 

*♦♦3.3  Volumetric  Scaling  Relations.  For  most  of  the  geometric 

shapes  which  best  describe  craters  (e.g. , paraboloid,  hyperboloid), 

2 

volume  is  proportional  to  r . It  is  frequently  assumed  that 

("°l)d  ’ ("V?°2)v  ’ (V”3)r  <M> 

where  n^  , ng  , and  n^  represent  the  scaling  exponents  for  depth, 
radius,  and  volume,  respectively.  Intuitively,  it  would  seem  that 
n 2 £ 1.0  , although  apparent  violations  have  frequently  been  noted. 
Development  of  scaling  relations  for  crater  volumes,  in  addition  to  the 
problems  outlined  above,  is  plagued  by  difficulties  in  obtaining  good 
volume  measurements.  Volumes  given  in  the  literature,  especially  those 
obtained  prior  to  I960,  are  oftentimes  based  on  one  or  two  radial  sur- 
veys and  are  therefore  not  as  accurate  as  those  developed  from  aerial 
stereophotography  or  from  a number  of  radial  surveys. 

4.3.4  Depth-o f-Burial  Scaling.  Cube-root  (mass  or  mass-gravity) 
scaling  ic  used  in  this  report  for  charge  DOB,  a common  approach.  Two 
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alternatives  are  also  found  in  other  studies: 

1.  If  the  appropriate  scaling  exponent  is  known  or  assumed  before* 
hand,  it  can  also  be  applied  to  COB  scaling. 

2.  An  iterative  approach  can  be  used  in  which  yield  is  held  con* 
stant  and  an  exponent  is  found  which  best  matches  both  crater  dimension 
and  DOB  (Reference  6). 

Generally,  all  three  methods  produce  satisfactory  data  fits;  some 
advantage  may  be  noted  in  the  last  method  where  data  are  grouped  to 
facilitate  its  use.  This  was  not  always  the  case  in  this  study,  how- 
ever. A primary  purpose  oJ  this  report  was  to  develop  and  illustrate 
dimension  scaling  relations;  hence  the  choice  of  cube-root  DOB  scaling. 
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Figure  U.ll  Composite  graph  for  apparent  crater  depth  for  1-ton 
TFT  spheres. 
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Figure  U.19  Composite  graph  for  apparent  crater  depth  for  1-kt  charges 


.on*  versus  charge  yield,  shoving  the  variation  of  scaling  ex- 


CHAPTER  5 


CRATER  EJECTA 

Although  not  a primary  consideration  in  this  report,  the  ejection 
mechanism  is  sufficiently  important  in  the  formation  of  a crates'  to 
warrant  explanation.  A complete  treatment  is  properly  the  subject  of 
a separate  study;  for  this,  the  reader  is  referred  to  the  Bibliography. 
The  following  paragraphs  are  Intended  to  provide  a general  appreciation 
of  this  phenomenon. 

5.1  COMPOSITION  AND  ORIGIN 

Crater  ejecta  consists  of  that  portion  of  the  soil  or  rock  debris 
thrown  beyond  the  boundaries  of  the  apparent  crater  by  an  explosion 
(Figure  5.1).  Together  with  the  fallback,  which  lies  between  the  true 
and  apparent  crater  boundaries,  it  comprises  all  material  completely 
dissociated  from  the  parent  medium  by  the  explosion.  It  may  represent 
a significant  hazard  at  considerable  distances  from  the  crater. 

In  a buried  explosion,  ejecta  mainly  originates  In  the  dome  of  earth 
material  which  is  forced  upward  by  the  expanding  gas  bubble.  As  the 
dome  disintegrates  and  gas  venting  occurs  (Figure  5*2),  discrete 
particles  of  ejected  material  enter  a trajectory  which  is , except  for 
very  small,  windborne  particles,  essentially  ballistic.  For  near- 
surface explosions,  where  the  fireball  obscures  the  cratered  region, 
the  mechanics  of  ejection  are  not  so  well  known.  Observations  of  near- 
surface  explosions  show  an  early,  fast-moving  corona  of  material 
ejected  from  a position  near  the  charge  and  at  a steep  angle  to  the 
ground  surface.  The  ejection  process  is,  however,  known  to  take  place 
over  a longer  period  of  time  and  to  include  lower  exit  angles. 

Material  fractured  by  the  ccetpressive  stress  wave  may  be  dislodged 
and  ejected  by  the  explosion  gases,  as  visualised  in  Figure  5.3.  How- 
ever, attempts  to  predict  ejecta  ranges  by  consideration  of  shock-front 
conditions  and  by  calculation  from  early  trajectory  parameters  have  been 
unsatisfactory.  It  appears  that  additional  experimental  observations 
will  be  necessary  for  this  purpose. 


-~c  r*}*  U- 


Theoretical  studies  indicate  that  in  large  near-surface  detonations, 
ejecta  particles  of  considerable  size  may  be  captured  by  the  thermal  up- 
drafts and  lofted  into  high,  nonballistic  trajectories.  Figure  5*^  il- 
lustrates origins  and  relative  ranges  for  general  cases  of  HE  charge 
geometries,  as  deterjjdned  from  field  observations. 


5.2  DESCRIPTIVE  PARAMETERS 

The  ejecta  field  is  divided  into  two  zones:  (l)  the  crater  lip 

(the  continuous  ejecta  surrounding  the  apparent  crater),  and  (2)  the 
discontinuous  ejecta,  comprising  the  discrete  natural  missiles  falling 
beyond  the  crater  lip. 

The  principal  parameters  used  to  describe  the  ejecta  are  the  aver- 
age lip  crest  height  (h)  (Figure  1.2);  the  radial  extent  of  the  crater 
lip  (r^)  from  GZ;  the  depth  of  deposition,  ejecta  mass  density,  and 
missile  size/dlstrlbution  defined  as  functions  of  radial  distance  from 
GZ;  and  the  maximum  missile  range  (r#).  The  principal  variables  which 
control  the  ejecta  parameters  are  the  shot  yield  and  geometry  and  the 
physical  nature  of  the  earth  medium. 

5.2.1  Crater  Lip.  The  amount  and  extent  of  the  continuously 
deposited  ejecta  in  the  crater  lip  are  determined  primarily  by  the  shot 
yield  and  geometry.  The  radial  extent  of  the  crater  lip  will  usually 
vary  from  about  2 to  4 apparent  crater  radii.  The  maximum  depth  of 
ejecta  in  the  lip  occurs  at  or  near  the  lip  crest,  and  its  height  above 
original  ground  can  be  estimated  as  about  one-fourth  to  one-third  the 
apparent  crater  depth  for  near-surface  burst*.  For  deeper  bursts,  the 
lip  height  is  usually  one- fifth  to  one- fourth  the  apparent  crater  depth. 
The  depth  of  ejecta  will  decrease  rapidly  in  an  exponential  fashion  as 
the  distance  from  GZ  Increases.  In  general,  the  volume  of  ejecta 
deposited  within  the  crater  lip  varies  from  about  Uo  percent  to  over 
90  percent  of  the  total,  the  latter  figure  representing  near-optimum 
DOB's  in  cohesive  media.  Deeper  bursts  may,  of  course,  result  in  a lip 
containing  all  of  the  ejected  material,  surrounding  an  apparent  crater 
of  insignificant  size  (Figure  4,22e).  Figure  5.5  illustrates  the 
fraction  of  total  deposition  with  range  from  GZ. 
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| ^ . 5.2.2  Discrete  E,1ecta  Field.  The  discontinuous  ejecta  beyond  the 

| crater  lip  is  usually  described  by  areal  mass  density  (e.g. , pounds  of 

■j  ejecta  per  square  foot)  and  a numerical  density  (number  of  missiles  per 

square  foot).  Both  parameters  decrease  exponentially  as  distance  from 

• 

GZ  increases,  and  a wide  circumferential  variation  is  usual,  which  is 
mainly  the  result  of  medium  inhomogeneities.  Figure  5.6  provides  a means 
of  predicting  ejecta  areal  density  for  a given  explosion.  Predictions 
obtained  from  these  curves  should  be  considered  only  as  first-order 
approximations , as  deviation  can  be  caused  by  variations  in  shot  geom- 
etry, earth  media  characteristics,  and  the  asymmetry  of  the  ejecta 


field  Itself. 

1 

Ejecta  missile  size  and  quantity,  both  in  the  crater  lip  and  in  ; 

the  discontinuous  portion  of  the  ejecta  field,  depend  primarily  on  the 
characteristics  of  the  earth  medium.  For  example,  the  cohesiveness  of 
a soil  with  a high  clay  content  often  results  in  missiles  of  substantial 
size,  while  a noncohesive  material  such  as  sand  will  produce  almost  no 
missiles  of  significant  size.  Explosions  in  glacial  tills  will  produce 
a large  n’rnber  of  long-range  missiles.  In  rock,  the  spacing  of  Joints 
is  a controlling  factor  in  determining  missile  size.  At  this  writing, 
means  of  predicting  natural  missile-size  distribution,  either  analyti- 
cally or  empirically,  are  considered  too  tentative  for  inclusion  in  a 
report  of  this  nature. 

For  buried  charges,  maximum  missile  range  is  approximately  pro- 
portional to  *1/6 

, and  is  shown  graphically  as  a function  of  DOB  in  I 

Figure  5«T  (from  Reference  7).  Surface  or  near-surface,  aboveground 

0 3 

bursts  produce  maximum  missile  ranges  which  more  nearly  scale  as  V . 

It  has  been  observed  that,  for  near-surface  HE  geometries,  the  periphery 

» 

of  the  ejecta  field  consists  of  predominantly  small  (<1  pound)  ! 

particles. 

Finally,  the  reader  is  cautioned  that  the  figures  of  this  chapter 
are  primarily  qualitative,  showing  trends  and  very  rough  values  for  -j 

ejecta  parameters.  Bo  attempt  should  be  made  to  extract  from  them  j! 

detailed  quantitative  information.  § 
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Figure  5-1  Throvout  of  ejecta  by  a low-yield  cratering 
exploaion  at  near-optimum  depth  of  burial. 


CHAPTER  6 


EFFECTS  OF  VARIATION  IN  SHOT  GEOMETRIES 
6.1  MULTIPLE-EXPLOSION  ARRAYS 


Weapons  or  munitions  may  be  detonated  in  close  proximity  to  one 
another  for  the  purpose  of  creating  military  obstacles  or  for  excava- 
tion purposes  (Figure  6.1).  When  the  detonations  are  simultaneous  and 
sufficiently  close  to  permit  interaction  betveen  charges,  the  shot  geom- 
etry is  herein  termed  "multiple  explosion."  A linear  array  of  this 
nature  is  termed  a "rov  shot,"  and  the  resulting  crater  is  referred  to 
as  a "rov  crater"  or  a "channel."  The  channel  may  also  result  from  tvo 
or  more  connecting  rov  shots.  The  rov-shot  geometry  has  been  the  object 
of  videspread  research,  due  to  its  numerous  possible  military  and  civil 
applications.  Multiple-explosion  arrays  may  also  be  tailored  to  other, 
less  common  purposes,  hovever,  and  may  be  nonlinear  and/or  detonated  in 
varying  sequences  . A portion  of  the  Bibliography  is  devoted  to  this 
general  cratering  application. 

6.1.1  Rov  Craters.  In  addition  to  the  cratering  variables  dis- 
cussed in  Chapters  3 aud  4,  rov  shots  are  dependent  upon  the  spacing 
betveen  charges  and  the  degree  of  simultaneity  of  detonation.  Sig- 
nificant departure  from  detonation  simultaneity  may  degrade  row-charge 
crater  dimensions.  Close  spacings  of  charges  (less  than  1.4  single- 
charge crater  radii)  enhance  both  crater  radius  and  depth,  as  compared 
to  single  craters  at  the  same  scaled  DOB.  Figure  6.2  shovs  the  en- 
hancement of  single -charge  crater  dimensions  obtained  by  reducing  the 
spacing  betveen  charges  in  a rev  charge  at  optimum  DOB.  As  the  charge 
spacing  decreases,  the  optimum  DOB  for  each  charge  must  be  increased 
by  the  enhancement  factor.  The  length  L of  a rov  crater  can  be 
found  from  the  equation 

L • s(N  - 1)  ♦ Zt%  (6.1) 

Where:  H ■ the  number  of  charges  in  the  rov 
s ■ charge  spacing 
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Spacing*  of  about  3rft  in  soil  and  up  to  l*r^  in  rock  may  produce 
satisfactory  linear  obstacles,  since  the  crater  lip  is  part  of  the 
obstacle. 

6.1.2  'Other  Multiple-Explosion  Arrays.  As  explained  in  Chapter  2, 
a variety  of  geometries  has  been  studied,  mostly  in  granular  soil,  for 
specific  applications  of  multiple-explosion  arrays.  These  have  included 
nonlinear  arrays,  adjacent  rov  shots  fired  simultaneously  and  in  varying 
sequences,  and  "multiple-pass'*  geometries,  wherein  row  charges  are  fired 
beneath  channels  created  by  preceding  row  shots.  The  applications  in- 
clude explosively  formed  earth  dams,  overburden  removal,  and  the  shaping 
of  row-shot  channels  to  certain  specifications.  Since  these  applica- 
tions are  quite  specialized,  no  attempt  will  be  made  in  this  report  to 
discuss  them  further.  Those  having  such  an  interest  should  consult 
the  Bibliography. 

6.2  CHARGE  STH4MIHG 

"Stemming”  refers  to  the  backfilling  of  material  in  the  charge- 
emplacement  hc’e.  Ideally,  charges  should  be  completely  stemmed  and 
tamped  to  a density  equal  to  that  of  the  parent  material  to  insure  that 
the  explosive  provides  the  maximum  in  cratering  performance.  There  are, 
however,  occasional  requirements,  mostly  military,  that  stemming  be  re- 
duced or  that  its  emplacement  and  removal  be  expedited.  Thus,  some 
attention  has  been  given  to  the  effects  on  crater  dimensions  of  various 
sizes  of  emplacement  holes  with  different  depths  of  stemming  and  with 
different  stemming  materials,  including  water.  Based  on  HE  experiments, 
it  has  been  concluded  that  stemming  does  little  to  increase  crater 
diameter,  which  is  generally  the  dimension  of  greatest  military  impor-' 
tance.  Stemming  of  about  one-half  the  emplacement-hole  depth  (50  percent 
stemming)  provides  most  of  the  crater  depth  which  would  be  expected  frem 
a fully  steamed  charge.  Water  appears  to  be  an  efficient  stemming  mate- 
rial. Figure  6.3  illustrates  the  small-scale  experimental  results,  as 
compiled  in  Reference  8. 
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6.3  UNDERWATER  CRATERING 

A limited  amount  of  experimentation  haa  been  devoted  to  cratering 
underwater.  With  one  exception  known  to  the  authors,  all  such  experi- 
ments within  the  time  frame  of  consideration  in  this  report  have  in- 
volved charges  resting  on  or  only  partially  buried  beneath  the  .earth- 
water  interface.  The  exception  is  Project  Tugboat  (by  EERL) , for  which 
both  small  and  large  HE  tests  of  buried  charges  have  been  conducted. 

The  goal  of  this  study  is  a light-draft  harbor  at  Kawaihae  on  the  island 
of  Hawaii. 

Reference  9,  which  is  included  in  the  Bibliography  on  this  subject, 
is  in  itself  a compendium  on  underwater  cratering  preceding  Tugboat. 

Data  from  a number  of  tests  are  compiled,  interpreted,  and  analyzed  in 
this  reference.  In  addition  to  single  charges,  limited  data  are  in- 
cluded on  row  and  "nail-driving"  detonations  (see  Section  6.6).  The  data 
from  the  a ingle -charge  experiments  are  included  in  Figure  6.4,  seeded  to 
the  1-ton  yield  level.  In  general,  both  row  and  nail-driving  experience 
parallels  that  for  land  craters,  with  allowances  made  for  differences  in 
single-charge  crater  shapes  and  sizes. 

Two  additional  5-ton  TNT  charges  fired  at  Mono  Lake,  California,  in 
1966  resulted  in  underwater  mounds  rather  than  craters.  Due  to  the 
anomalous,  unexplained  results,  these  shots  were  not  included  in 
Reference  9.  Also,  a small  quantity  of  old,  formerly  classified  data 
has  recently  been  located,  and  this  may  be  included  in  future  analyses. 

6.4  DEEPLY  BURIED  EXPLOSIONS 

The  cratering  effects  of  explosive  charges  buried  deeper  than  opti- 
mum DOB  are  illustrated  in  Appendix  A and  in  Figures  4.1  through  4.19. 

A continued  increase  in  DOB  will  result  in  a smaller  apparent  crater, 
although  true  crater  dimensions  will  continue  to  increase.  At  some  DOB, 
the  apparent  crater  will  cease  to  be  evident,  or  may  (especially  in  co- 
hesive material)  result  in  a mound  due  to  heaving  and/or  bulking  action. 
Essentially,  this  is  the  containment  DOB.  Below  this  depth,  a eamouflet, 
or  underground  cavity,  is  formed.  The  true  crater  no  longer  Intersects 
the  surface , but  is  coincident  with  the  eamouflet,  with  concentric  zones 
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of  deformation.  These  conditions  are  illustrated  in  Figure  4.22. 

Approximate  containment  depths  may  he  extrapolated  from  the  figures 

of  Chapter  4 in  terms  of  DOB/W*/3  . The  depth  of  containment  vas 

1/3 

reported  in  Reference  2 as  shout  3.5  ft/lb  . Cavity  radii  are  ap- 
proximately 1.2  ft/lb^3  for  HE  charges,  or  about  47  ft/kt^^3  for  HE,  as 
taken  from  reports  on  a number  of  nuclear  experiments.  For  HE,  an 
empirical  equation  vhich  considers  DOB  is 

re  * ~‘'DOB~ — (6.2)  (Reference  10) 
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More  accurate  calculations  for  NE,  vhich  take  into  account  the  effects 
of  medium  properties  and  DOB,  provide  the  equation 

y1/3 

r ■ C — *—  (6.3)  (Reference  11) 

C (pz)° 

Where:  C * an  observed  cavity  proportionality  constant 

a <■  an  exponent  based  on  the  adiabatic  (gas)  expansion 
coefficient 

Other  terms  are  as  previously  defined  in  the  text  or  Figure  1.2.  The 
explosion  gases  are,  of  course,  formed  by  vaporization  of  the  medium. 
Figure  6.5  shows  values  of  a as  a function  of  vat er  content,  and 
assumed  values  of  a , p , and  C are  shown  below  for  several  rock 
media. 


Medium 

a 

P 

C 

gm/cm3 

Granite 

0.324 

2.7 

103 

Tuff 

0.292 

1.9 

97 

Alluvium 

0.296 

1.9 

89 

Salt 

0.311 

2.3 

96 

Dolomite 

0.329 

2.3 

89 

To  use  Equation  6.3,  enter  W In  kilotons  and  Z in  meters;  cavity 

radius  r will  be  found  in  meters, 
c 

In  desert  alluvium,  subsidence  craters  have  been  experienced  for 
large  yields  fired  below  containment  depths.  This  action,  which  may 
take  place  over  several  days,  results  from  the  subsidence  of  overlying 
soil  ir.to  the  cavity  formed  by  the  explosion.  Characteristically, 
these  craters  are  wide  and  shallow.  In  view  of  the  uncertainties  as- 
sociated with  their  formation,  no  prediction  techniques  appear 
applicable. 

6.5  BOMB/SHELL  CRATERS 

As  indicated  in  Chapter  2,  cratering  by  conventional  bombs  has  long 
been  of  interest  to  the  military.  Recently,  limited  cratering  tests 
have  been  conducted  with  artillery  and  mortar  shells  in  connection  with 
the  design  of  protective  structures.  It  does  not  appear  sufficient  to 
predict  bomb/shell  crater  dimensions  from  those  formed  by  bare  charges 
equivalent  in  yield  to  the  bomb  or  shell  filler  material,  even  though 
it  may  be  argued  that  explosive  energy  lost  in  rupturing  the  case  is 
regained  in  kinetic  energy  of  the  fragments.  The  explosive  filler  ma- 
terial is  usually  cast  in  a cylindrical  configuration,  and  detonation 
is  usually  initiated  in  the  nose  or  tail  of  the  warhead;  these  consider- 
ations also  fail  to  answer  completely  the  questions  raised  over  differ- 
ences in  craters. 

The  limited  study  which  has  been  done  on  this  subject  indicates 
that  deeply  buried  bombs  produce  craters  smaller  in  radius  bat  larger 
in  depth  than  do  comparable  bare  charges.  There  do  net  seem  to  be  suffi- 
cient data  on  shallowly  buried  bomba  to  permit  conclusions.  Experience 
with  mortar  and  artillery  shells,  which  contain  a proportionately 
heavier  casing,  indicates  that  near-surface  (fuze  quick)  bursts  form 
craters  which  are  somewhat  larger  than  those  of  bare  charges.  Unfortu- 
nately, a complete  resolution  of  the  problem  is  hampered  by  lack  of 
definitive  test  data.  Figure  6.6  shows  general  bcob-data  trends,  and 
may  be  compared  with  spherical-charge  curves  in  Appendix  B. 
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6.6  SUCCESSIVE  SHOTS  ON  A VERTICAL  AXIS 

Some  interest  haa  been  shown  in  the  craters  produced  by  a aeries 
of  surface  or  near-surface  explosions,  each  detonation  (after  the  first) 
occurring  in  the  center  of  the  crater  left  by  the  preceding  shot.  This 
technique,  known  as  nail  driving,  may  represent  a means  of  attacking 
deeply  buried,  hardened  structures.  The  results  of  several  experiments 
conducted  in  granular  and  clayey  soils  and  in  competent  rock  are  shown 
by  envelopes  of  data  in  Figures  6.7  and  6.8.  Although  there  is  consid- 
erable scatter  in  the  data,  the  envslopes  approximate  the  expected  in- 
creases in  crater  dimensions , and  also  show  the  general  differences  to 
be  expected  in  soil  and  rock. 
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Figure  6.3  Increase  in  HS  crater  dimensions  as  functions  of 
stemming  and  DOB. 


SCALED  WATER  DEPTH.  FT/Lb’'* 

Figure  6.k  Apparent  dimensions  of  underwater  craters  from  1-ton  charges  fired  at  the  earth- 
water  interface  in  a variety  of  fine-grained  materials.  Anomalous  appearance  of  curves  in 
very  shallow  water  (left  side  of  graph)  is  attributed  to  washback  effect,  in  which  ejected 
material  is  carried  back  into  crater. 


CHAPTER  7 


ENVIRONMENTAL  INFLUENCES 

The  graphs  in  Appendix  B and  Chapter  l*  are  based  primarily  on  tests 
conducted  in  relatively  homogeneous , isotropic  media.  In  many  situa- 
tions , however,  it  may  be  necessary  to  detonate  explosives  in  a medium 
containing  some  environmental  anomaly,  such  as  a water  table  at  a shal- 
low depth,  a layering  of  one  type  of  soil  over  another,  parallel  planes 
of  distinct  Jointing  in  rock,  or  a steeply  sloping  ground  surface.  All 
of  these  factors  can  influence  the  formation  of  a crater  and,  in  some 
cases,  radically  change  its  size  or  characteristics. 

7.1  SLOPING  TOPOGRAPHY 

Terrain  slopes  of  about  5 degrees  or  greater  will  affect  crater 
formation  for  a surface  explosion,  the  venting  process  for  a buried  ex- 
plosion, and  ejecta  distribution  in  any  case.  For  gentle  slopes,  the 
total  volumetric  effects  will  be  about  the  same  as  for  craters  on  level 
ground,  but  the  resulting  crater  will  be  asymmetrical,  wider  up-slope 
and  with  a larger  lip  down-slope  (Figure  7.1).  For  the  field  of  discrete 
ejecta  particles,  greater  maximum  ranges  will  occur  down-slope,  assuming 
that  the  wind  is  not  a significant  factor. 

Limited  small-scale  cratering  experiments  have  been  conducted  in 
moist,  sandy  soil  and  in  desert  alluvium  on  slopes  ranging  from 
•*0  degrees  to  vertical  wedges,  the  latter  representing  the  extreme  in 
sloping  topography.  For  charges  buried  on  severe  but  nonvertical  slopes, 
with  DOB  measured  from  the  sloping  surface  and  with  the  vertical  depth 
of  overburden  being  greater  than  containment  depth,  crater  dimensions 
decrease  with  increasing  slope.  Optimum  vertical  DOB  and  optimum  dis- 
tance from  the  free,  vertical  face  of  a wedge  appear  roughly  the  same, 
and  perhaps  larger  by  about  one-third  than  optimum  DOB  on  level  terrain. 
For  this  geometry,  ejecta  distribution  is  preponderantly  directed  toward 
the  free  face,  with  about  three-fourths  of  the  total  ejecta  mass  falling 
in  this  direction  when  DOB  is  optimum  for  crater  volume.  The  disparity 
in  ejecta  distribution  increases  with  further  increase  in  DOB. 
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7.2  -LAYERED  SYSTEMS 
# 

7.2.1  Water  Tables.  Based  upon  HE  experiments , a subsurface  water 
table  In  a soil  medium  will  begin  to  influence  the  size  and  shape  of  the 
crater  when  its  depth  below  the  surface  is  equal  to  or  less  than  three- 
fourths  the  predicted  apparent  crater  radius.  Its  effect  is  to  flatten 
and  widen  the  crater.  As  the  water  table  depth  decreases,  its  effect 
becomes  more  evident;  for  a water  table  depth  of  one- fourth  to  one- 
fifth  the  original  predicted  crater  radius,  the  final  radius  may  be  as 
much  as  50  percent  greater  than  and  the  depth  as  little  as  one-third 
that  of  the  original  predicted  value. 

7.2.2  Bedrock.  The  influence  of  a bedrock  layer  below  a soil 
medium  is  similar  to  that  of  a water  table,  though  somewhat  leas  pro- 
nounced. For  HE  explosions  at  the  surface,  the  bedrock  layer  may  in- 
crease the  crater  radius  slightly  (5-10  percent),  and  may  decrease  the 
final  depth  by  as  much  as  one-third  when  the  overburden  layer  is  as 
shallow  as  one-fourth  the  predicted  apparent  crater  radius. 

7.2.3  Rock  Bedding /Jointing.  For  low-yield  NE  and  high-yield  HE 
explosions  at  or  very  near  the  surface,  the  bedding  or  jointing  planes  in 
rock  can  Influence  the  shape  of  the  crater  produced,  as  well  as  the  di- 
rection of  the  ejection  process.  The  formation  of  the  crater  will  tend 
to  follow  the  direction  of  the  predominant  joints,  thus  increasing  the 
crater  radius  by  as  much  as  one-third  in  the  direction  parallel  to  the 
joints,  or  decreasing  it  by  as  much  ao  one-third  normal  to  the  joints. 

The  magnitude  of  the  crater  depth  is  usually  not  affected  significantly, 
but  the  deepest  point  may  be  shifted  to  one  side  of  the  crater.  As  yield 
or  DOB  is  Increased,  the  influence  of  rock  Jointing  is  reduced. 

The  dip  of  bedding  planes  will  influence  energy  propagation, 
causing  the  maximum  crater  depth  to  be  offset  in  the  down-dip  direction. 
Little  overall  effect  is  noted  in  regard  to  crater  radius,  but  differ- 
ences in  ejection  angles  cause  the  maximum  lip  height  and  ejecta  radius 
to  occur  down-dip. 
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7.3  SNOW  AND  ICE 

Measured  craters  in  snow  or  ice  are  a rarity,  and  for  this  reason 
are  presented  separate  from  the  more  general  data  of  Appendix  B and 
Chapter  1*.  A few  craters  have  been  recorded  for  surface  HE  explosions 
in  snow/ice;  these  are  larger  than  craters  in  soil,  and  are  characteris- 
tically wide  and  flat.  Figure  7.2  shows  trends  in  crater  size  and  shape 
for  a surface-hurst  geometry  in  this  medium. 
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(Courtesy  of  Lawrence  Radiation  Laboratory 
and  0.  S.  Atomic  Energy  Commission) 


CHAPTER  8 


o 1 

SUMMARY  AND  RECOMMENDATIONS  I 

. i 

\ 

3.1  SUMMARY  » 

The  foregoing  chapters  have  presented  a brief  history  of  cratering 
research  and  a compilation  of  HE  and  NE  crater  data.  An  attempt  has 
been  made  to  sort  these  data  in  such  a fashion  as  to  identify  the  aig-  i 

nificant  variables  affecting  crater  size,  and  an  analysis  has  been  con- 

\ 

ducted  to  provide  the  tools  necessary  for  crater  predictions.  Varia- 
tions in  crater  shapes  and  sizes  due  to  departures  from  the  usual  test 
geometries  or  to  unusual  environmental  conditions  have  been  discussed. 

Phenomena  associated  with  debris  ejection  from  the  crater  have  also  been 

i 

I 

examined  briefly.  The  result  is  a series  of  graphs  from  which  crater 
parameters  may  be  read  directly  (Appendix  B)  or  from  which  data  trends 

can  be  identified  (Chapters  U and  5)»  and  to  which  certain  Judgment  • 

1 

factors  can  be  applied  for  practical  usage.  j 


8.2  RECOMMENDATIONS 


The  report  has  been  prepared  in  a manner  to  facilitate  the  addition 
of  crater  data  by  the  user.  It  is  recommended  that  a formal  updating  be 
accomplished  on  a periodic  basis — say,  biennially  or  briennially — to 
insure  maximum  use  of  new  data  as  it  becomes  available.  With  the  data 
retrieval  and  computerized  plotting  system  which  has  been  established, 
this  could  be  done  without  great  effort.  For  the  more  specialized 
crater  applications  (e.g.,  multiple-charge  arrays,  cratering  on  slopes), 
such  a periodic  updating  would  help  fill  serious  gaps  in  the  existing 
data.  To  aid  in  updating,  users  are  urged  to  submit  suggestions  and 
corrections  which  they  feel  are  appropriate. 

A crater/ejecta  study  should  be  a part  of  every  experimental  plan, 
and  especially  so  for  tests  involving  NE  " large  HE  yields.  In  the 
past,  where  this  has  been  overlooked  or  neglected,  the  result  has  been 
the  loss  of  data  needed  at  some  later  time.  In  general,  the  area  in 
which  crater  data  are  n»st  needed  is  that  of  large  yields  in  clays  and 

92 


BSS3iS5BEOBS!a52Baa5SaS 


sassssas 


V silts  and  mixtures  thereof  which  profcabl y represent  the  major  portion  of 

the  earth's  soils.  More  information  is  needed  on  cratering  in  layered 
systems,  also,  to  insure  proper  application  of  available  data  to  practi- 
cal problems.  True  crater  dimensions  and  zones  of  subsurface  deformation, 
often  omitted  from  explosion  test  research,  should  be  measured  wherever 
possible.  These  data  are  important  not  only  in  the  prediction  of  damage 
to  underground  structures,  but  also  in  the  formulation  of  volumetric 
and  mass-balance  equations,  as  well  as  in  formulating  expressions  for 
basic  cratering  mechanisms.  Accurate  volumetric  data,  particularly  of 
true  craters,  are  valuable  for  normalizing  certain  other  energy  input 
phenomena,  e.g.  ground  motion,  cratering  efficiency,  etc.  Finally, 
every  available  opportunity  to  study  the  phenomenon  of  crater  ejecta 
should  be  exploited,  since  it  is  the  damage  mechanism  in  cratering  which  ■ 

is  potentially  the  most  far  reaching.  Considerable  research  is  needed 
to  accurately  quantify  the  parameters  discussed  in  Chapter  5 • 

} 
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TABULATION  OF  CRATER  DATA 


This  appendix  contains  the  data  tabulation  (Tables  A.l  through 
A.lU)  introduced  in  Section  3.3.  All  but  Table  A.lU  are  reproductions 
of  computer  printouts,  the  program  for  which  is  contained  in  Appendix  D. 
Explanatory  material,  also  included  in  the  program,  is  reproduced  below. 

IDENTIFICATION  OF  CODES  AND  VARIABLES 


STC  Operation  Stagecoach 

RRV  Railroad  Vulnerability  Program 

SES  Cratering  by  Ground  Burst  at  Suffield  Experimental  Station 

SCH  Project  Schooner 

PB  Project  Pre-Buggy 

AV  Air  Vent  Cratering  Series 

MS  C Mine  Shaft  - Calibration  Cratering  Series 

MS  I-MO  Mine  Shaft  I - Mine  Ore 

MS  I -MU  Mine  Shaft  I - Mine  Under 

MS  II-MR  Mine  Shaft  II  - Mineral  Rock 

BKBD  Operation  Buckboard 

SC  II  Sandia  Corporation  Cratering  Series  II 

SC  I Sandia  Corporation  Cratering  Series  I 

MTCE  Multiple  Threat  Cratering  Experiment 

DB  Project  Danny  Boy 


DP  Operation  Distant  Plain 

TPOTS  Teapot  Ess 

PAL  Palanquin 

FT  Flat  Top 

ST AMS  Simulation  Testa  of  Artillery  and  Mortar  Shell  Explosions 

J,JS,JU  Operation  Jangle 

CSM  Colorado  School  of  Mines,  Underground  Explosion  Test  Program 

UETP  Underground  Explosion  Test  Program 

JSSC  Cratering  in  Sand  from  Spherical  Charges 
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EUE  Effects  of  Underground  Explosions 

FCBT  Fort  Churchill  Blast  Tests 

MOLE  Small  Explosion  Tests,  Project  Mole 

CESBC  Cratering  Effects  of  Surface  and  Buried  HE  Charges 

CDS  Cratering  in  Dry  Sand 

ESRIC  Effects  of  a Soil-Rock  Interface  on  Cratering 

SEU  Effects  of  Stemming  on  Underground  Explosions 

SE  Stemming  Effects  for  Certain  HE  Charges 

PCE  Energy  Partitioning  for  Partially  Confined  Explosions 

ICSBA  Cratering  Tests  in  Basalt  - Inter-Oceanic  Canal  Study 

ICSCC  Cratering  Tests  in  Cucaracha  Culehra  - Inter-Oceanic  Canal 

Study 

ICSGS  Cratering  Tests  in  Gatun  Sandstone  - Inter-Oceanic  Canal 

Study 

ICSMM  Cratering  Tests  in  Marine  Muck  - Inter-Oceanic  Canal  Study 

ICSRC  Cratering  Tests  in  Residual  Clay  - Inter-Oceanic  Canal  Study 

FICS  Investigation  of  Charge  Shape  at  Ft.  Churchill 

SRI  Crater  Study,  Operation  Castle,  Stanford  Research  Institute 

ANC  Ammonium  Nitrate  Cratering 

HUS  Russian  Nuclear  Event- 

PG  Project  Pre-Gondola 

NEP  Project  Neptune 

TRIN  Project  Trinidad 

JJ  WES  Stemming  Series 

PF  Operation  Prairie  Flat 

SEDAN  Operation  Sedan 

SL  Sandia  Laboratories  Series 

ZUL  Project  Zulu 

DIAL  Operation  Dial  Pack 
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MEDIUM  CCDE 


Concrete  or  Grout  0 

Hard  Rock-Granite,  Basalt,  Limestone,  Etc.  . 1 

Soft  Rock-Sandstone,  Etc 2 

Very  Soft  Rock  or  Very  Hard  Soil  - Shale, 

Tuff,  Frozen  Soil 3 

Clay It 

Silty  Clay 6 

Loess  and  Lacustrine  Silt 5 

Silt,  Sand,  and  Clay 7 

Silty  Sand  and  Desert  Alluvium 8 

Sand 

MOISTURE  CODE 

Coiqpletely  Dry 0 

Very  Dry 1 

Dry 

Slightly  Moist 3 

Moist 4 

Very  Moist 5 

Slightly  Vet 6 

Vet 7 

Very  Vet 8 

Saturated  9 

EXPLOSIVE  CODE 

Nuclear 0 

TNT  . . 

C-U,  C-3 2 

PentoJLite 3 

Ammonium  Nitrate  - Slurry  or  Grains  4 
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Nitromethane 5 

Dynamite 6 

Amatol 7 

C-3  and  Tetrytal*  8 


CRATER  SHAPE  CODE 
0-Unknown 
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(Continued) 
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CRATER  SHAPE  CODE  (Concluded) 


6-PAN 


7-CONE 


9-MOUND 
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TABUS  A.l4  CRATER  DATA  FOR  HEMISPHERICAL  CHARGES 


Event 

Yield 

(TNT) 

Med- 

iun 

Code 

Mois- 

ture 

Code 

Apparent  Crater 

Radius 

Depth 

Up 

Haight 

Volune 

pounds 

feet 

feet 

feet 

ft3 

Sandy  Silty  Clay: 

SES*  (5-*hot  series, 

512 

7 

3 

6.2 

4.8 

1.0 

mm 

Fall  1558) 

551 

7 

3 

6.1 

3.8 

0.5 

m m 

523 

7 

3 

5.3 

3.0 

0.3 

mm 

521 

7 

3 

5.8 

3-3 

0.4 

mm 

520 

7 

3 

6.2 

3.4 

0.2 

mm 

SE3  (Fall  1958) 

600 

7 

3 

6.6 

3.7 

0.5 

mm 

SES  (Oct  1959) 

10,000 

7 

3 

20.0 

12.5 

mm 

mm 

8ES  (Aug  I960) 

40,000 

7 

3 

36.7 

15.1 

— 

m m 

SES  (Aug  1961) 

200,000 

7 

3 

69.9 

20.7 

— 

mm 

SES  (Jul  1963) 

40,000 

7 

3 

:e.i 

21.3 

— 

mm 

SES  (Aug  1963) 

40,000 

7 

3 

40.7 

19.7 

— 

mm 

8ES  (Aug  1963) 

10,000 

7 

3 

19.7 

14.1 

mm 

mm 

SES  (Sep  1963) 

10,000 

7 

3 

20.1 

14.1 

mm 

mm 

Snowball 

1,000,000 

7 

3 

139.8 

13.8 

— 

580,000 

Distant  Plain  4 

100,000 

7 

3 

40.7 

16.4 

— 

— 

Dry  Sand: 

White  Tribe  1-1 

H,f6C 

9 

1 

lM 

6.6 

2.1 

mm 

White  Tribe  1-2 

11,560 

9 

1 

20.0 

10.3 

2.0 

— 

White  Tribe  1-3 

11,560 

9 

1 

16.0 

6.2 

2.1 

mm 

White  Tribe  II-l 

11,560 

9 

1 

18.3 

8.6 

1.8 

mm 

White  Tribe  11-2 

11,560 

9 

1 

17.9 

8.3 

2-3 

mm 

White  Tribe  n-3 

11,560 

9 

1 

18.0 

9.6 

2.6 

mm 

White  Tribe  III-l 

11,560 

9 

1 

17.3 

6.1 

1.5 

— 

White  Tribe  III-2 

11,560 

9 

1 

15.5 

7.5 

1.9 

a*  m 

White  Tribe  III-3 

11,560 

9 

1 

19.1 

10.0 

1.6 

mm 

* Conducted  at  Suffleld  Experimental  Station  (now  Defence  Research  Estab- 
lishment, Suffleld),  Alberta,  Canada. 


APPENDIX  B 


GRAPHICAL  PRESENTATION  OF  CRATER  DATA 

The  tabular  data  of  Appendix  A are  contained  in  graphical  form  in. 
Figures  B.l  through  B.93  of  this  appendix.  Essentially,  these  are 
least-squares  fits  of  the  data,  as  explained  in  Section  3.3.  .'Where  in- 
sufficient data  exist  to  support  a curve  (actually,  a strai.y'.t-line  fit 
in  all  but  Figure  B.93),  the  available  points  are  ahovn.  Whs re  no  data 
exist  for  one  or  tvo  plots  in  a figure  (apparent  dimensions,  true  dimen- 
sions, or  volumes},  the  blank  graph(s)  is  (are)  Included  for  future 
entries.  This  permits  the  use  of  a consistent  format,  vith  all  graphs 
in  a single  figure  printed  and  bound  for  ease  in  reading. 

Where  available,  .classified  data  have  been  considered  in  prepara- 
tion of  the  graphs,  although  not  ahovn  here  or  elsewhere  in  this  re- 
port. Existence  of  such  data  are  noted  on  the  appropriate  figures. 

The  category  numbers  used  in  the  graphs  and  in  the  tables  in 
Appendix  A refer  to  the  height-of -burst  categories  into  which  the  data 
are  divided  (as  explained  in  Section  3*3). 
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Figure  B.l  Dimensions  of  craters  in  basalt  and  granite  for 
0.20  < Z < 0.50  ft/lb1/3,  Category  2 (sheet  1 of  2). 
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k TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 
figure  B.3  Distensions  of  craters  in  'basalt  and  granite  for 
-0.05  < Z < 0.05  ft/lb^i  Category  k (sheet  1 of  2). 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHAROE  YIELD 

Figure  B.4  Dimensions  of  craters  in  basalt  and  granite  fur 
-0.20  < 2 < -0,05  ft/lb1^,  Category  5 (sheet  1 of  2). 
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Figure  B.4  (sheet  2 of  2) 
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Figure  B.5  Dimensions  of  craters  in  basalt  and  granite  for 
-0.50  < Z < -0.20  ft/lb1^3,  Category  6 (sheet  1 of  2). 
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Figure  B.6  Diaenaiona  of  cratera  la  % an  alt  and  granite  for 
-0.90  < Z < -0.50  ft/lb1^  Category  7 (ataeet  1 of  2). 
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Figure  B.7  Dimoslons  of  craters  la  basalt  and  granite  for 
-1.10  < S < -0.90  ft/lb1^3,  Category  8 (sheet  1 of  2). 
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Figure  B.8  Distensions  of  craters  in  basalt  and  granite  for 
-2.00  < Z < -1.10  ffc/lb1^3,  Category  9 (sheet  1 of  2). 
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Figure  B.9  Dimensions  of  craters  in  basalt  and  granite  for 
Z < -2.00  ft/lb1^,  Category  10  (sheet  1 of  2). 
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Z < -2.00  Category  3.0  (ai«efc  .1  of  8). 
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Figure  B.16  Diaensioas  of  craters  In  ah  ale , tuff*  sod  frozen  growid 
for  0.05  i Z < 0.20  ft/lfc^,  Category  3 (sheet  1 of  ?). 
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c.  a*»are:<t  ano  true  crater  volumes  versus  charge  yield 
Figure  B.16  (sheet  2 of  2). 
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Figure  B.17  Dimension*  of  crater*  in  *hale,  tuff,  and  frozen  ground 
• for  -0.05  < Z < 0.05  ft /It1/3,  Category  4- (sheet  I' of  2). 
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rigure  B.19  >'lnenaions  of  craters  in  shale,  tuff,  and  frozen  ground 
for  -0.50  < 3 < ‘-0.20  ft/lh1^,  Category  6 (sheet  1 of  2). 
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Figure  8.21  Dimensions  of  orators  in  shale*  tuff,  and  frozen  ground 
for  -1.10  < Z < -0.90  Category  8 (sheet  1 of  2) 
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figure  B.22  Dimensions  of  craters  in  shale,  tuff,  and  frozen  ground 
for  -2.00  < Z < -1.10  ft/lb1/3,  Category  9 (sheet  1 of  2). 
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Figure  B.24  Dimensions  of  craters  in  dry  clay  for 
0.50  <.  Z ft /lb1/3,  Category  1 (sheet  1 of  2). 
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Figure  B.25  Dimensiono  of  craters  in  dry  cloy  for 
0,20  < Z < 0.50  ft/lb1//3,  Category  2 (sheet  1 of  2). 
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Figure  B.2 6 Dimensions  of  craters  in  dry  clay  for 
0,05  < Z < 0.20  ft /l b1^,  Category  3 (sheet  1 of  2). 
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Figure  B.27  Dinensions  of  craters  in  dry  clay  for 
-0.05  <,Z  < 0.05  ft/lb^3,  Category  U (sheet  1 of  2). 
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Figure  B.31  Dlcenniono  of  craters  in  dry  clay  for 
-1.10  <.  Z < -0.90  ft/lb1/3,  Category  8 (sheet  1 of  2), 
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Figure  B.32  Dimensions  of  craters  in  dry  clay  for 
-2.00  < Z < -1.10  fw/lb1^3,  Category  9 (sheet  1 of  2). 
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Figure  B.33  Dimensions  of  crater*  in  dry  clay  for 
Z < -2.00  ft/lb1/3,  Category  10  (sheet  1 of  2). 
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Figure  B.'lU  Dimensions  of  craters  in  moist  clay  for 
-0.05  <,  Z < 0.05  ft/lb1^,  Category  4 (sheet  1 of  2). 
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Figure  B.36  Dimensions  of  craters  in  moist  clcy  fo- 
-C.90  <.  Z < -0.50  ft/lb1^,  Category  7 (sheet  1 of  <?), 
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Figure  B.37  Dimensions  of  craters  in  moist  clay  for 
-1.10  < Z < -0.90  ft/lb1/3,  Category  8 (sheet  1 of  2). 


Figure  b. 37  (sheet  2 of  2). 
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Figure  B.l'O  Dimensions  of  craters  in  wet  clay  for 
0.05  <.  Z < 0.20  ft/lb1^,  Category  3 (sheet  1 of  2). 
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Figure  B.U6  Dimensions  of  craters  in  wet  clay  for 
Z < -2.00  ft/lb1^3,  Category  10  (sheet  1 of  2). 
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Figure  B.U6  (sheet  2 of  2). 
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Figure  B.l*7  Dimensions  of  craters  in  moist  loens  and  moist  lacustrine 
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b.  TRUE  CRATER  OlkT-NSIONI  VERSUS  CHARGE  YIELD 
Figure  B.U8  Dimensions  of  craters  in  moist  loess  and  moist  lacustrin< 
silt  for  -0.05  < Z < 0.05  ft/lb1^3,  Category  U (sheet  1 of  2) 
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Figure  B.l*9  Dimensions  of  craters  in  moist  loess  end  moist  lacustrine 
silt  for  -0.20  < Z < -0.05  ft/lb1^3,  Category  5 (sheet  1 of  2). 
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Figure  B.50  Dimensions  of  craters  in  moist  loess  and  moist  lacustrine 
silt  for  -0.50  < Z < -0.20  ft/lb1^,  Category  6 (sheet  1 of  2). 
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h.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  3.51  Dimensions  of  craters  in  moist  loess  and  moist  lacustrine 
silt  for  -0.90  < Z < -0.50  ft/lb1^3,  Category  7 (sheet  1 of  2). 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  B.52  Dimensions  of  craters  in  moist  loess  ani  moist  lacustrine 
silt  for  -1.10  £ Z < -0.90  ft/lb1^3,  Category  8 (sheet  1 of  2). 
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Figure  B.55  Diaeniions  of  craters  in  ooist  silty  clay  for 
0.50  < Z ft/lb^^,  Category  1 (sheet  1 of  2). 
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Figure  B.5&  Dimensiona  of  craters  in  moist  silty  clay  for 
-0.50  <_  Z < -0.20  ft/lto1^,  Category  6 (sheet  1 of  2). 
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lx  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 


Figure  B.6l  Dimensions  of  craters  In  dry-to-aoist  sandy  silty  clay 
for -0.05  <,2  < 0.05  ft/ib1^,  Category  4 (sheet  1 of  2). 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  B.62  Dioenoions  of  craters  in  desert  alluvium  for 
0.05  < Z < 0.20  ft/lb1/3,  Category  3 (sheet  1 of  2). 
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Figure  B.63  (sheet  2 of  2) 
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a.  APPARENT  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 
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k TAU3  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  B.6U  Dimensions  of  craters  in  desert  alluvium  for 
-0.20  < 2 < -0.05  ffc/11)1^3,  Category  5 (sheet  1 of  2). 
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h.  7RU6  CRATER  DIMENSION*  VERSUS  CHARfl*  YIELD 


Figure  B.66  Dimensions  of  craters  in  desert  alluvium  for 
-1.10  < Z < -0.90  ft/ro1/3,  Category  3 (sheet  1 of  2). 
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e.  APPARENT  AND  TRUE  C:V,TER  VOLUMES  VERSUS  CHARGE  YIELD 

Figure  B.69  (sheet  2 of  2). 
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e.  ARRAREftT  AND  TRUE  CRATER  VOLUMES  VERSUS  CHARGE  YIELD 

Figure  B.70  (sheet  2 of  2), 
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k TRUE  CRATER  DIMEXO'lOW*  VERSUS  CHARGE  YIELD  • 

Figure  B.71  Dimeasio.  a of  craters  in  moist  sanciy  silt  for 
-0.05  < Z < 0.05  ft/lV1^3,  Category  4 (sheet  1 of  2). 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 


Figure  B.T2  Dimensions  of  craters  in  aoi3t  sandy  silt  for 
~0.90  < Z < -0.50  ft/ib"  , Category  7 (sheet  1 of  2). 
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V'tAuY  CRAnYoiMlNsTowf VERSUS  CHARG8  YIELD 


Figure  B.73  Dinenaionu  of  craters  in  soiat  sandy  silt  for 
-1.10  < Z < -0.90  ft/lb1/3,  CVegory  9 (sheet  1 of  2). 


LEGEND 

APPARENT  CRATER  RAOIU3 
APPARENT  CRATER  DEPTH 
APPARENT  CRATER  UP  HT 


CRATE*  DIMENSIONS.  FEET  CRATER  DIMENSIONS.  FECT 


CHflRCC  YIElD.  POUNDS 

•.  APPARENT  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 


_ 1 
" 5 

o 

X 

" I 

" 5 

o 

x 

~ I 
- S 


• 5 
o 


- 5 
o 

X 

~ 1 
7 S 


NO  OATA  AVAILABLE  {FOR  FUTURE  USE) 


LEOENO 

0 

TRUE  CRATER  RADIUS 

1 

TRUE  CRATER  DEPTH 

_!  1 1 nil* 

LilMIh!  LUlllUi  1 1 

i 1 IIIIS 

5 l 

5 1 SI  5 1 

5 

I SI  SI  SI  5 1 St 
tXlO  * I 1X10  1 I (XIO  * I (XIO  * I 'X10  4 ) !XI0  * I (XIO  * ) !XJ0  ’ I (XIQ  * I !X10  * I 

CHARGE  TIElO.  '0UN05 

b.  TRUE  CRATER  DIMENSIONS  VERCUS  CHARGE  YIELO 


•Figuie  B.76  Dimensions  of  craters  in  dry-to-xoint  sauid  for 
0.20  ' Z < 0.50  ft./lb1^3,  Category  2 (sheet  1 of  2). 
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t,  AFPARENT  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELO 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  B.7H  Dimensions  of  craters  in  dry-to-noist  sand  for 
-0.05  <.  Z < 0.05  ft/lt^3,  Category  4 (sheet  1 of  2). 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 


Figure  B.T9  Dimensions  of  craters  in  dry-to-moist  sand  for 
-0*20  < 1 < -0*05  ft/lb1^3,  Category  5 (sheet  1 of  2). 
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Figure  B.79  (sheet  2 of  2) 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELO 

Figure  B.80  Dimensions  of  craters  in  dry-to-moist  sand  for 
-0.50  <^Z  < -0.20  ft/lh^-^,  Category  6 (sheet  1 of  2). 
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a.  ARfARENT  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 
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k TRUE  CRATEn  DIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  B.8l  Dimensions  of  craters  in  dry-to-moist  send  for 
-1.10  <Z<  -0.90  ft/lb1^3,  Category  8 (sheet  l of  2). 
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i 


1 


CHARGE  YIELD.  P0UN03 

k TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 


Figure  B.82  Dimensions  of  craters  in  dry-to-moist  sand  for 
-2.00  <.  Z < -1.10  ft/lb1^3!  Category  9 (sheet  1 of  2). 
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k.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  B.03  Dimensions  of  craters  in  vet  sand  for  0.50  < Z ft/lb^ 
Category  1 (sheet  1 of  2). 
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k TRUE  CRATER  DIMENSIONS  VERSUS  CHAROE  YIELD 

Figure  B.8U  Dimensions  of  craters  in  vet  sand  for 
0.20  < Z < 0.50  ft/lb1^,  Category  2 (sheet  1 of  2). 
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c.  AFRaRENT  and  true  crater  volumes  versus  CHARGE  YIELD 


Figure  B.8U  (sheet  2 of  2) 
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4.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  VIILD 
Figure  B.85  Dimensions  of  craters  In  vet  sand  for 
0.05  <.  Z < 0.20  rt/lb1^3.  Category  3 (stoet  1 of  2). 
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b.  TRUE  CRATER  OIMENSIONS  VERSUS  CHARGE  YIELD 

Figure  B.66  Dimensions  of  craters  A a vet  sand  for 

-0.05  * Z < 0.05  ft/lb1^3,  Cavegory  4 (sheet  1 of  2). 
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b.  TRUE  CRATER  DIMENSION*  VERSUS  CHARGE  YIELD 


Figure  B.87  Dimensions  of  craters  in  wet  sand  for 
-0.20  < Z < -0.05  ft/lb1/3,  Category  5 (sheet  1 of  2). 
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b.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARGE  YIELD  , 

Figure  B.83  Dimensions  of  craters  in  vet  sand  for 
-0.50  <,  Z < -0.20  ft/lb1^3,  Category  6 (sheet  1 of  2). 
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k.  TRUE  CRATER  DIMENSIONS  VERSUS  CHARCE  YIELD 

Figure  B.89  Dimensions  of  craters  in  vet  sar.d  for 
-0.90  < Z < -0.50  ft/lb1^3,  Category  7 (sheet  1 of  2).  . 
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Figure  B.93  Apparent  crater  radius  and  depth  versus  charge  yield  for 
hemispherical  changes  in  sandy  silty  clay.  Both  curves  are  based  upon 
15  data  points.  The  depth  curve  is  approximated;  for  yields  up  to 

0 Us 

about  10,000  pounds,  its  equation  is  0.20  W . 
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This  appendix  contains  a listing,  generally  in  the  order  of  intro- 
duction in  the  text,  of  reference  material  which  was  drawn  upon  most 
heavily  to  obtain  crater  data  or  to  prepare  the  synopsis  of  cratering 
research  in  Chapter  2 or  the  discussions  in  Chapters  5,  6,  and  7*  It 
is  divided  into  general  fields  of  cratering  subjects,  and  it  will  be 
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APPENDIX  D 


COMPUTER  PROGRAM  FOR  CRATER  DATA 

This  appendix  contains  an  annotated  listing  (Table  D.l)  of  the  com- 
puter program  developed  to  sort,  analyze,  and  plot  the  crater  data  in 
Appendixes  A and  B.  The  program  is  shown  in  this  form  in  the  hope  that 
it  may  be  more  understandable  to  the  layman,  vhile  at  the  same  time  pro- 
viding sufficient  information  from  which  a similar  program  could  be 
constructed  if  desired. 

Note  that  the  program  is  essentially  in  two  parts — a main  program 
and  a plotting  subroutine.  The  computer  language  is  FORTRAN^  IV;  the 
program  was  run  on  a GE-bOO  Series  computer  at  the  WES.  Plotting  was 
accomplished  on  a CalComp  Plotter. 


1 


Formula  translation. 
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TABLE  D.x  CCMPWER  PROOW  TOR  CRA TSR 


BEST  AVAILABLE  COPY 


1 

2 

1 

l 

5 

6 
1 

0 

* 

11 

11 

‘I? 

13 

1< 
i5 
.1* 
17 
1* 
l* 
20 
2t 
2? 
2S 
2* 
2’ 
20 
-2  7 
2« 
2* 
30 


a 

I 

i 

« 

3 

I 


31 

32 

33 

34 

35 
3» 
37 
3* 
3« 

40 

41 

42 

43 

44 

45 
4* 
47 
4* 
4* 


NOLml 

•coiK  

BCD0U7 

0JHENS10N  l®‘|^‘?®®>ev(2oo).MaB«200».**<200).OA<200>.MLt20#) 

COMMON  7S<7).TL<7>,E7l7>.1»'75> 
r*U  EcOTSUBUT  IlMOO®-3) 

REMIND  3 ,, 

fi*LL  ELOT (0 . • *30 . « *3 ) 

417  me*mns  effects  ° 

l|VtStON//”X‘»««*7£R  DATA  COMPILATION) 

1 fOBIlAT<t2*6.4l<.  H) *- 

tr(H.oT.o>  oo  T°  5i 

[MN.L>.01  CO  TO  52  . 

CO  TO  5 3 

* 

161  S«%.  <wO«0SnJMJ*l.I2> 

41>  f 0«H*T«lHl»25)l»12A47j! ,uv, 

40J  ^SItOH  .M;.3MO.H.?X.3HTOT.37.,HC06rr.7R.3B£^.  **.mi.7*.2HU 
1,*X.2HX7.10X.2HY2J 
M0*0 

pO  102  J*1 *10 

«0*0  ..  - — -• 

rjT*o. 

r«o. 

no  225  1*1 >7 

M(U*0 
SUNYIt  >*0  • 

sxrai«o.  - ■ 

St/HXTU  J *5* 

SX2IU«n. 

22>  CONTINUF 

»«»«» «.■»’ » ■">»’ 

IfTZU T!cE.<r.».*W'^?7^'  |0,2>  oo  TO  202 

j.EO.32  oo  TO  282 

00  T0  "* 


(03Bttnu«d) 
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BEST  AVAILABLE  COPY 

TABU!  D>1  (COTTINUED) 


Ixpla*  Card  Stateaent 

nation  Ho. 


>0 

>1 

>3 

54 

>9 

54 

17 


IF  1 2 1 1 ) . L T , >'0'.  05 . A^O  .'2 1 1 ) • GE '.  *9 . 2 • AUD . J « EG  • 9 ) GO  TO  202 


IT 

jf 
0 ' 
1 
2 
3 

FT " 
9 

*- 

7 

• 

9 

9“ 

1 

2 

3 

4 

79 


IF(Z(J).LT,-0.9.AND.Z(I ) ,GE .-1 .1 . AND.*1, 
!H2<I>.lT.-l.'i;AHDrZ<n.e£.-2.o;ANO.J. 
1FIZ(1>.IT.-2.0.AND.J.E0.10)  GO  TO  202 

GO  TO  201  ' 

202  NO*NO*l 

rrTRxrrnioTTTznTTn 

H( 1 ) *H( 1 )*1 

su"mj.30Myfn'*A'W5Tiarrn 

SXY(1> t$XY(l)*AlOG(CY( 1 5 ••ALOOIHAI 1 5 > 

■ 9 U ttU  1 Ti  S Uh  X ITT«ALfl6'i  ?YTT>T 

9X2(1 )>9X2(l)»AL0G(CY(t > >*AU0G(CY( l ) > 

~2innuTrn-;M;5.7-G8'  to1  212 

H(2)*h(2)«1 

■ ■ • suhy  1 2 * iiuKr't  2 r*iuoffTwrrn 

SXY (2) *SXY(2>*AL0G( CY 1 1 > ) *ALOO( DA <1)5 

' SUHX72J«SUHXl2)*ALlJT!TirifnTJ 

SX2(2)aSX2|2)«AL0G:CY(I )>*AL00(CY(I >) 
~j\i  mmiD.EO.O.)  GO  TO  2l3 


EO.I 

£0.0 


) TO 

i J 

to; 

t 

fTfl 

n 

to; 

t 

) 00 
) GO 


TO  202 
TO  202 


H(3)aH(3)«l 

suhy  i3) » sohy  c ir*AXOffVHnrrr 


nr 

77 

7» 

70 


SXYl3)*SXYI3)*AL00«CY(|)>aAU00<H|.(!)) 

' SUHXI3 ) >SUHX(3 )*»AL0GTffY IT) )" 

SX2|3)*SX2(3)*aL0G(CY ( t ) )*AL0G(CY(  t ) ) 

-JI3'inRTUlTE0.0-.rfij  TO  214 

H(*)*H(4>»1 

SUMYU)tSUHYr4)T*tOGTirrnT) 

$XY(4)tsXY(4).AL0G(CY(t )>*ALOG(RT(I>) 

• - 9 U H X (4  r « S U MXT4T*  A tff  47CVTI D 

SX2(4)«SX2(4)»AL0G(CY(t ))*AL0G(Cy(I)> 

714'  ir  CffTUT^0.T7I  GO-TO  2 ITS 

HI9)*H(9)*1  

’ SUHY  1 3 ) »$UHY|3 ) ♦ ALOOT&tTTT) 
9XY(9)«9XY(9)»AL0G(CY( I ) )»AU00(0T( I ) ) 
SUHXl5)*SUHX(9)iAi;00fCY'fI)>  — - - 
9X2(9 » »SX2(5>»AL00(CY ( 1 > >»ALOQ(CY( 1 ) > 
219  IMVA(  i >,to.O.)  GO  TO  211 
H(4)*H($)*1 

" ~suHf ( jtcuhy vatid 

SXY(4)>SXY(4>*A(.0G(CY(t))*AL0Q(VA(t)) 

soHx(*)«5onx<4>»AuOG(i;Y<Trr 

9X2(4) *9X2(4) *AL0a(CY( t )>*ALOG(CY(t>) 
-7n-Trro7-mT£a:o.)  go  m 211 

HI7)*H(7)*1 

' 9oHYi7)»suHYTrrvAi:«nmrrT) 

9XV(7)«9XY(7)»AL0GICY(  I ) )*A|>OOIVT(  I ) ) 

— swxm»s\mm»i\.mzrnrr 

9X2(7) •9X2(7)«ALOOICY(t  MaALOOICY(t)) 

*717’  COHTINU* 

lr(HUt). £0.0.0)  Ht(  1 )*0.001 

Y3I)I0)*A100I0  TTTCTTTT»7fl 

y3(N0)*y3(N0)*0.79 


ir* 

10 

16 

H 

12 

13 


09 

~16~ 

»7 

-^or- 


TOT* 
SI 


(Oontlnuad) 
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Hoi«h  UMtjytrM  Fit  C«lcul«tloca 


TABU!  0.1  (COBTimjED) 


Expla-  Ctrd 
flttlon  Ho. 


® available  copy 

Btutetunt 


“I 


*104 

105 


TIT 

107 
"100 
100 
~T10 
111 

“TIT" 
113 
“ITT 
115 
"Tl« 
117 

~ur~ 
u* 
Tao 
121 
122 
123 

T rr 

125 
120 
127 

m 

12* 

~mr- 

131 

132 

133 

~ur- 

135 
~I3*~ 
137 
130 
13* 
-1*0  - 
1*1 
~T«2" 
1*3 
"I**' 
1*5 
- I*« 
1*7 
TOO  — 


ir<vT<TKt0“O'rWT7Tff7WI- 

.T7<M)»»iooio»*TUn*i.o 

TTTNOr*Y7'('NIJJ«0 ,73 

irtOT(II.EO.O.O)  DTI | )»0.001 

T5mo)»*LooioTOTiTrr»7nr 

Y5  *N0>  *r5(N0 1«0 , 75 

irc«n  rv.ro-.  a.  aT-BTm»8;wr  - 

Y*(H0)Mt  OGlOlRTlt  ItO.fl 

— rrfwur»T*<MBr«0".75 


if<y*<i).ea.o.o>  vmi)»o,ooi 

Y6<NO)«iiosiO(viTrrRrro — - 

_ Y«(NO>»Y*(NO)*0.75 

i r « o*  1 1 > . rovoT  ■DTnurrirar 

Y2(NO)**L0010ID*|Fm»2.0 

YTTWrrrvTrTOTTJ-.TJ 

„ „ . . 1*<0AH  MjO.0.0)  (UUl'O.OOl 
Yl  fHO ) •>L3clO  (TfiTTTf  •TTiT — ~ 
YHNO)»Yl(H0>»0.75 

inno)«iBEH<n 

XlNOJ.ALOOlOlCYll)) 

•«nfo)«Tnnjrsj775 

IFIXlHOl.OT.n  fiXlNO) 

ir(x<NO).LT.rmTmifrKor- 

BEMT«r-rjT 

201  Iff l .LT.N1  GO  TO  103 

IflNQ.EQ.Ol  00  TO  10J 

T*nnrr37T- 


K*ExFioir/o.75) 

DO  210  K«1.7  •• 


ir<H|K}.LT.2)  GO  TO  40* 

{fOEHT.iT.-o.TirnBar'.o- 

EN»H(H) 


C*EN*SX2<K> 
D»3UHX(K  JiSUHYTiO" 
E»SUHXIHJ«SU«XIH) 

rr««-o--  

6«C*E 


irtir;cT;u.uuuui j mxi»o,o 

tflQ.LT. 0.00001)  GO  TO  *10 

"EXIKWT70 


*10  H»<SUHY<H).|x(K1«JWHXIIO>/EH 

COtH)*tXFf  HI 

Y5(K>»C0tK)»q,Q)««EnK) 


1** 

-T50- 

151 

152  — 

153 

~13* 

155 

“13* 

157 


GO  TO  *05 

*o*  Ystx)*o.3or 

YLtKKO.QOl 

cotn)«r. — 

EX(X>»0. 

“*03'CU*t IWUE 


H0*H0»1 

-jr  (HOvorrftr'GOTO'Tor 

•J03  COHTJHUE 


(Oootlautd) 
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Drat  First  Set  at  Axel 


tuts  o.i  (coittihued) 


BEST  AVAILABLE  COPY 


BxpU- 

Cud 

Ststaent 

Mtlon 

Do. 

SGJBOuTTnE'  BrATER 
REAL  LEG.IOENT 


COMMON  yS<7),yl<7),EX<7>.  |S(J5) 

0 1«EN5I0W"TL#C  TTJ  mirmTlEffT  2 r.TTJENTCS ) 
yLBUI)*  BhCRATER  0 

YCBnE^'lHTSfiNSTYJS 

YLiLtJ)*  8ms,  rEET 


CALL  LGAXISlO. ,0, « YLIL, *24 , ,1,  T.4.1, .01.1.2,1) 


CALL  PLOT (7, 9, 4,9. 2) 

' CALC  PlOTTT.SToVm 

CALL  LGAMS(«..0..*LBL.-2«..l.  L?  lMtl..t.2,i> 

“CALL"  PLOT  { 4779 V3«  3) " 

CALL  SLOT (4. . 1.4,2) 

“sill  piorrr,  .T:’4,2J 

CALL  PL0T(7. ,0.3,2) 

~ CALL-  Pton  47V  07TV21 

10ENT(l»t»MAPPARENT 

'10ENTC21MMCRATER 

|DEnT<3>»«hRA0IUS 


CALL  SYMBOL (9.2, 1.3. 0.1. LEG. 0, ,6) 

CALL  SYMBOL  f *7  A'.T.T.  OTl'.TOEKTi  0 . . 2« ) 

I DENT (2) >8m  CRATER 

njEHTT3T«BHBEPTH 

CALL  SYmB0l<4.«,0.B,0.1.10ENT.0..24) 


vmri 

10ERT(3»«BHL1P  mt 

"CALL  SYMBOL!  47*7 0757 a7ITTDFNT;07.2»1 

CALL  POINT (4, 2, 1.19.1,. 04) 

c*ivPoiirr(47r.-o7By^77o*T 

CALL  PO|NT|4,2.0.99,3..04) 


«*«  ) -41  ny 

IP(Yl ( I > ,LT . 0. > 00  TO  311 

irritTTr.T!T.-sr-oo'Tir3n 

CALL  P01NT(XM),Y1(I).1..«4) 

Sll'TTtYtTTl'.'LT.'OVT'OO  TO*  312 

lM|B(i>.ar.0>  GO  TO  312 


tel 


fr.:  «. 

vf=,  5^'  > .. 

A"  1 4’*' 

/ CfX  '* 


tabu  o.i  (coirrmiEo) 


IxpU*  Card 
utlon  Us. 


8t*t*aent 


C S *9 
o 5 41 
4 I *2 

g f-u- 

c a *4 

; *9 
■3  44 

*7 


ir  ■! 


I l \i 

■if« 
. ? ♦» 

*3 

* | 100 
* 2‘ioi' 


tl04 
"107 
101 


220  CONTINUE  - -.  .. 

lr(TS(2).LT.0.0)  60  TO  221 

•|rtgltt2)TamT5T~'i!0~I0  22T 
irCEKCEJ.EO.O.J  00  To  221 
CALL  OLOTCO.,T$(2),J) 

CALL  7L0T<r,YL<2>.2> 

221  CONTINUE 

|MYS<3).tT.0.0)  qo  TQ  222 


ir<E*(3).E0.0.)  00  TO  222 
CALL  0L0T(0.',TSI3>;3)' 

CALL  NLOm*.YLC3>,2> 

222  CONTINUE  

CALL  PLOTlO, ,0. >>3) 

~C<Lt;i6ATmff.tB;,VLBL,«24..i.  T,l 

CALL  *LO?C0. ,4.9,3) 

CALL  0L0T{7.y.4.f,2) 

CALL  OLOTC7.9.0..2) 

CALL  LOA* ISC 0 0 . . *lSL ;=24T.m ,~ 
CALL  Nc.0T(4. ,0.3.3) 

CIU'T»Ltrm77TTTT5 

CALL  0L0TC7., 1.0,2) 

’ ' CALL  *L0TC7.,07J,2J 

CALL  OLOTCO. ,0.3,2) 

-..CALL  SYH80L(5.2.1.3.D.I7tECr0.,*) 
WENT ( 1 )*0hTRUE  CNa 

1 Bg^ryTipTEiriAin 

|DENT(3)»0MU$ 

- - CALL  *TMBOLC4;i.0;O76nTrDINr."0'..2 

IDENT (2>*0HTEft  OEOT 
|DEnT(3)*Shh  ‘ 

CALL  STNOOLCO.O.O.S.O.l, I0ENT.0. ,2 
Wlf  ?OTNTI4,J,0,95,2,.0*) 

CALL  00|NT(4.3,0.99,3..00) 

‘ 00  302  t.iiNa 

ir<Y4C|).LT,0,»  00  TO  313 

irC|«C|).0T.0')"OO  TO  313 

CALL  00|NTCX(|),T4(!),2..00) 

"7iTTnT9iTr;irr.o;>  so  ia  302 

irciocD.or.o)  oo  to  302 

“CALC  Faprnx < IT;t5 «Ti . J.T«n 

302  CONTINUE 

■ tr<TSl4T.-LT;r.OT10  t5‘223 

inEX(4).0T.1.9)  00  TO  223 

rncrroriEo.oTrco  to  a;; 

CALL  0L0T(U.,TSI4>.3> 

call  olo r c r . Tl'crr, n 

223  CONTINUE 

irtrjC5»;LTr«:or'6iTTy-jm 

|f«EX|»).OT.l.»)  00  TO  22* 

irrrxTsr-.'iaiii.)  so  to  2» 

CALL  0L0T(0.,TS(9),3) 

- ~-CAtf9lorcriVI-{9T7rj 

220  CONTINUE 
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BET  AVAILABLE  COPY 


TAS1X  D.l  (CONCLUDED) 


Expla-  Cant 
nation  Mo. 


Stateacnt 


J 


TUUjjatwsr'cuiTr 

yLBL(4)i  8M  FEET 


mi'nsrroTH.r^yj 

CAU.  LGAX|S<0.*0..YL8L*»32».1»  T 1 • 2.1J 

CALL  FLOT(Or;A.“97TT 

CALL  PL0T<7.5.4.3,2> 

CALL"  PL0TrTr57ir77TJ 


CALL  LOAXIStO* >0. iXLBL.*24# . 1> 

TXLC>UJTir.-,0’.3;3J 

CALL  PLOT ( 4 , , 1 , A > 2 ) 

CALC  PL0TLT;7r7U72T 

CALL  PL0Tc7.,0.3.2> 

call  pLOTi4:.a;r,^Y 

CALL  SYM80L<5.2.1.3.0.1>LE3.0..«> 

T0Pnrr»-iT»rtFPCTPn 

tOENT (2> *8'l  CRATER 

iDEMTa)-.vHvaLuHf 

CALL  STHBOL(4.«.0.«.0.1.!OENT.0..24l 

I0EnT<  it  *8hTR‘uECAa 

IOENT|2)alHTER  volu 

nTERTT3TJlftKJ 

CALL  JYNBOL44.6.0.5.0 .1. 10ENT.O. >24 1 

CALL'  P 0 1 N Y<  4 - 2.0 . W7  2 r.7T*1 

CALL  POiNT<4.2,0.35.3..0B) 

DO  303  m.Nfl  

irmm.LT.o.)  eo  to  sir 


T«  10 1 0 » 1 ■ » li  2»  1 > 


CALL  P01NT{X(I).T«(J).2..0»I 
irfrTdi.cr.o'.rTJirTO'SOT — 
iruit n.oT.oi  ao  to  303 


call  pointer < n rmn.x.;3tr 

CONTINUE 


i 

. ». 

. -m 

5 H2 

rrTT34ff)  .TT  .B’.O  L 5U  'T0  '224 — 

!F<£X(*).GT.1.5>  00  TO  224 

a 

< 

9 143 

S 144 
i 147  - - 

3 14* 

TP4PH*) .E070T) "GITTO  224 

CALL  PLOT(0.,YS<«).3> 

CALL  PtOTTr.TLT*J«‘2T 
224  CONTINUE 

V T *T“ 
A 14B 

IPl  V$<7) *LT*0*flT  AO  10  229 
!HiX(?)«GT.l»5)  GO  TO  225 

t**  — 

150 

CALL  PLOTIo! » Y$< 7)#S> 

• 

X5l  "" 
152 

CAUL  PL0TT7  « YU  7 ) #21 
225  CONTIHUC 

l 

r “ tjj 

$ i»« 

C*LU  PLOT 1 17*  » -SO*  $ *37 
CALL  PLOTCO*  #3* ••3) 

1 

llM' 
0 »» 

return  - - 
END 
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Cratering  programs  and  data  resulting  from  numerous  single-charge  explosion  tests 
are  Siam  cited  and  compiled  in  tabular  font.  Analyses  are  performed  on  these  data 
to  provide  means  of  predicting  basic  cratering  parameter*.  Prediction  equations 
are  developed  from  least-squares , straight-line  plots,  and  are  presented  in  graphical 
form.  Means  of  updating  these  tabulations  and  analyses  on  a regular  basis  by  auto- 
matic data  processing  are  discussed.  Data  I'l  grouped  so  as  to  account  for  the 
factors  which  primarily  affect  crater  sise  and  shape:  yield,  burst  geometry,  and 
cratered  medium.  The  Influence  of  other  conditions,  such  as  soil  moisture,  layered 
media,  etc.,  is  also  considered,  tophssis  is  on  single-charge,  dry-land  experiments, 
which  best  permit  Isolation  of  the  factors  contributing  to  the  basic  parameters. 
However,  effects  of  environmental  influences,  unusual  charge  geometries,  and  other 
factors  significantly  affscting  craters  are  also  briefly  considered.  Similarly, 
baste  ejects  phenomena  are  included.  Trends  in  crater  dimensions  are  shown  by  metns 
of  graphs  normalised  to  charge  sixes  cosoeaturablt  to  large  chemical  and  small 
nuclear  yields.  Sealing  a a a prediction  tool  la  discussed. 
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